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Summary 


The equivalent widths of Hy have been determined for 137 Southern 
B-stars (tabulated in Section 5), using material obtained from the Radcliffe 
Observatory in Pretoria (on the average 3 spectra per star). Most of the 
spectrophotometric reductions were carried out at the Dominion Astrophysical 
Observatory in Victoria, B.C. Detailed analysis, described in Sections 3 and 4, 
showed that there were no systematic differences, either between different 
observers, or between the present southern work and the previous investiga- 
tions of northern stars by Petrie at Victoria. ‘The internal error derived 
for the present determinations of Hy intensities is 7 per cent per star, 
equivalent to o™-23 in the absolute magnitudes. ‘To convert these absolute 
magnitudes into stellar distances, a uniform system of star-magnitudes and 
colour-indices has been established for the southern hemisphere (Section 7). 
The preceding Section 6 discusses the luminosity-spectrum array. 

Section 8 provides satisfactory tests of the newly-determined spectro- 
scopic distances, in particular with the help of mean proper-motion parallaxes 
and, furthermore, by using mean distances derived from galactic-rotation 
effects on the radial velocities. In the latter case, it is shown (Section 8.4.3), 
that we may not be justified in using the south “‘ calcium-cloud distances ”’ (i.e. 
‘* effective half-distances ’’, based on the K-line) since, differing from the 
north, lines of sight in the south mostly run through cloud-free inter-arm 
regions. 

Section g relates the new material to the structure of the Milky Way, and 
makes use of recent work on H I-areas, clusters, associations, H 11-regions 
and various proposed Milky-Way arm-patterns. A general colour-excess 
diagram, illustrating the dependency on galactic longitude, is derived in 
Section 9.3. Figs. 9, 10 and 11 represent the B-stars and other Milky-Way 
constituents in their relation to different galactic regions out to about 4, 1°5, 
and 7 kiloparsecs, respectively. Particular attention is drawn to the 
occurrence of a conspicuous group of reddened B-stars which are clearly 
situated between hydrogen arms, lying especially outside the Sagittarius arm. 
Various significant discrepancies between location of stars, neutral hydrogen, 
and rotational models of the Galaxy are pointed out. Evolutionary aspects are 
indicated. 





1. Introduction.—‘‘ It has long been recognized that one of the main difficulties 
in galactic studies is the uncertainty in the distance of the B-stars. Direct 
Measurements are too few and uncertain . . . we are therefore compelled to develop 
the method of spectroscopic absolute magnitudes and thus to assign distances 
to the individual B-stars...’’. There can hardly be a better introduction to this 
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paper than this remark made by Petrie early in his work on the determination of 
absolute magnitudes of B-stars by measurement of the total intensity of the 
hydrogen lines. Petrie’s papers (1953, 1956a, 1956b, 1956c, 1958a, 1958b, 
1958c, 1960) give bibliographical notes on work in this field by other authors. 
Special attention may also be drawn to the recent publications by Kopylov 
(1958a, 1958b) and Sinnerstad (1961). 

The fundamental basis of this type of work is the existence of a close correlation 
between the total strength of Hy and the absolute magnitude of the star in whose 
spectrum it occurs—a relation which has been calibrated by standard stars of 
known luminosities. By combining absolute magnitudes inferred in this way 
with other data (apparent brightness, colour, proper motion, radial velocity) 
we can derive the individual distances of these stars and check at least the statis- 
tical accuracy of the results. All work of this kind has so far been restricted to 
B-stars accessible to northern observatories. However, an extension to the 
southern sky has been made possible by the use of spectroscopic material from the 
Radcliffe Observatory, Pretoria, obtained by Feast, Thackeray, and Wesselink 
(1955, 1957, 1958), primarily for measuring B-star velocities. Dr Thackeray has 
very kindly put this material at the writer’s disposal. 

2. The southern B-star spectra.—The Radcliffe Observatory’s B-star programme 
was specially designed to be complementary to the northern B-star work that has 
been carried out for many years at the Dominion Astrophysical Observatory, 
Victoria, B.C., first by J. S. Plaskett, who was later joined by J. A. Pearce, and 
then more recently by R. M. Petrie. 

Although the Radcliffe and Victoria spectrographs are of quite different design 
—in the arrangements in which they have been used for most of this work—they 
give very nearly the same dispersions, so that methods successfully developed at 
Victoria could be applied with very little modification directly to the Pretoria 
material. ‘The Pretoria spectrograms were obtained with two cameras, the 
b-camera giving 29 A/mm at Hy, and the c-camera giving 49 A/mm; the spectra 
belong to the Pretoria ‘‘OBD”- and “‘ OBC ”’-series which contain mainly 
fainter stars. The plates used were Kodak I1aO or 103aO emulsions; on each 
one were two calibration exposures, made with a tube-sensitometer, supplying 
12 spots which were arranged so as to minimize the effects of emulsion irregularities. 
The spectrograms were intended primarily for radial velocity work, not for spectro- 
photometry, and this must be remembered when judging the errors of the line- 
intensities. Doubtless somewhat greater accuracy would have been obtained, 
had it been possible to introduce all the special precautions of accurate photo- 
graphic photometry. 

Between 1 and g spectrograms were available for each of 137 southern B-stars, 
with an average of three spectrograms per star. ‘The stars ranged down to 10™-2 
(visual). Fig. 1 shows their galactic distribution; the data for the southern stars 
come from Table III (Section 5), and for the northern stars are those referred to in 
Section 9.1. The galactic coordinates used here are /' and b! (Ohlsson 1932); see 
also Blaauw, Gum, Pawsey, and Westerhout (1960). 


3. Spectrophotometry and equivalent widths of Hy. 

3.1. The measurement of the plates, as far as the main programme was con- 
cerned, was carried out with the Cambridge microphotometer, using tracings at a 
magnification 100: 1, and covering about 100A centred on Hy. The calibration 
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spots were recorded at1o:1. ‘The tracings usually included the A4388A He -line, 
although normally this line lay well outside the wings of Hy. Later work showed 
that occasionally Hy has very extensive wings, making it necessary to record a 
greater length of the spectrum—up to 200A—to ensure a better determination of 
the continuous background, and a number of plates were re-photometered on this 
account. 

3.2. Several other lines, mainly O 11, may occur in appreciable strength within 
the Hy-wings (see, e.g., Struve, 1931; and Butler and Seddon, 1958). In 
allowing for these lines, and also in determining the line wings, the same general 
procedure was used as had been followed in the original Victoria work, in the hope 
that in their error-characteristics the new southern measures would be closely 
similar to the earlier northern measures. 

3-3. In outline, the procedure was to measure the equivalent width of Hy, 
to read off the corresponding absolute magnitude from the calibration-curve 
(Petrie 1953, p. 264) and finally to add a correction for spectral type (/.c., p. 263). 

Preliminary Cambridge results for a few stars showed that satisfactory accuracy 
was being achieved with relatively few spectrograms. For instance, 8 spectro- 
grams of HD 149438 gave a mean equivalent width 3:62A with a standard error 
+015, or 4 per cent. Three spectrograms of HD 93128 gave 1-51 + 0°16, i.e. 
an error of 11 per cent. ‘These and other samples indicated that a mean of only 
3 plates would usually give an equivalent width with a standard error less than 
1opercent. The general aim, therefore, was to use not less than three spectrograms 
per star. 

3.4. However, before the work had gone very far, Dr Petrie offered very helpful 
assistance. He kindly arranged that the writer should spend four months at the 
Dominion Astrophysical Observatory, so that the very efficient semi-automatic 
intensity-reduction machine there could be used (Beals 1944). The visit to 
Victoria also permitted much closer consultation with Dr Petrie, and a more 
thorough intercomparison of all Victoria and Pretoria results available up to date. 

3.5. In examining the errors of the Hy-intensities we first note that, if N is the 
number of stars for which m plates were measured, the distribution n(V) of the 
observations falls into the following m groups: 1(8), 2(33), 3(58), 4(32), 5(3) 
6(1), 8(1), 9(1). The values of m are given in Column 11 of Table III. The 
average number of plates per star is 3-02. There are 129 stars with n>2, 
i.e. altogether 7 groups, which can contribute to the error-examination. 

The standard errors were first computed from internal agreement, using the 


formula 
W/ zr 
AW= © n(n—1)’ (1) 


where the 5-values are individual deviations from the mean of nm measures. ‘The 
AW-values (in units of o-o1 A) are given in Column 13 of Table III. 

Two diagrams were then drawn, one of which was a plot of the absolute values 
AW against W, the other of the percentage errors AW (per cent) against W. These 
showed that AW tended to vary with W, and that for some purposes it is better to 
deal with the percentage error in W rather than with the absolute error. If, then, 
A; is the percentage-deviation of the equivalent-width determination (derived 
from the ith plate) from the mean (based on z plates), we find the r.m.s. deviation 
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for each star by 


2= 1 S(ay. (2) 
—T ga 
The mean value of o? per plate, averaged over all the N stars in one group, each 
having » plates per star, is therefore 


2— 2 
=a 2 ],: (3) 
The values of o,? were then derived for each group using the above formulae. 
Finally, the mean [o?] over all 7 groups was formed by weighting each one according 
to Nn. 

From the above-mentioned 129 stars we thus obtained [o?] = 66755/406, i.e., 
[o]= +12°8 per cent; so that the average star, with 3-02 plates, has an error of 

+7°4 per cent in its Hy-intensity. 

We may convert this Hy error into the equivalent error in absolute magnitude, 
using Petrie’s equation M=a.10°", where a= — 10-96 and b= —0-145. It will 
be shown in Section 4 that we are entitled to apply this equation to the Pretoria 
material, although the formula was originally derived for the northern stars only. 

3-6. Were there a constant 7-4 per cent error in the Hy intensity, the errors in M 
would range between about o™-2 and o™-3, according to the value of M. We have, 
however, computed errors in more detail for Table III, to try to give some indica- 
tion of the relative reliabilities of the M-values. By differentiating Petrie’s 
formula, we obtain 

AM .log,,e=b.AW. M, 
and 
AM/M = —0°335.AW. (4) 


Table III gives the AW-values computed from the interagreement of different 
spectrograms; its Column 15 contains the corresponding A from this formula, 
expressed in units of o™-o1. The values of M used are, however, not quite 
those given in the preceding Column 14 of Table III, but those which were derived 
before the spectral-type correction was applied (see Section 3.3). The mean 
value of all the AM in Table III is o™-23. 

Petrie (1953, p. 265; 1960) has given the standard error of the absolute 
magnitudes in his work as o™-3. He has also remarked that the scatter of indivi- 
dual points about the calibration-curve iso™-5. This scatter may arise from errors 
in the assumed parallaxes of the standard stars, or from an intrinsic scatter in the 
(Hy, M)-relation. From whichever source this standard scatter arises, the 
internal consistency of both the Victoria and the Pretoria work appears to be quite 
satisfactory in view of our present knowledge of B-star luminosities. 


4. The correspondence between Pretoria and Victoria luminosities. 

4.1. It is of the greatest importance that there should be as little systematic 
difference as possible between Pretoria and Victoria magnitudes, and that the 
error-characteristics of the one should closely resemble those of the other, so that 
the results may be used for the study of B-stars equally all round the Milky Way. 

4.2. The first important question is whether there are systematic differences 
between different observers when dealing with the same material, differences such 
as may arise from individual judgments as to the position of the continuum, etc. 
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Dr Petrie kindly repeated 22 of the author’s measurements and found that for 
B-stars his measured Hy equivalent-widths were on the average about 10 per cent 
smaller than the author’s, and that for O-stars the differences were somewhat 
larger, possibly due to variation in the interpretation of the blended lines lying 
within Hy. A 10 per cent difference would correspond to about o™-3 in absolute 
magnitude. Later, however, the author was able to increase the overlapping 
material between the two observers to 46 spectrograms of 10 stars, and the result 


was then: mean difference AW (Beer-Petrie) = — 0-02 A; the difference was still 
less if each star was weighted according to the number of spectrograms involved. 
The values for individual stars are shown in Fig. 2. It might be mentioned that 


there appears to be a slight tendency for the AW (Beer-Petrie) to increase with 
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Fic. 2.—Individual differences in measured equivalent widths of the same stars as obtained by 
different observers; following the No. of the star (Table III) is the number of plates (in brackets). 


photographic density in the continuous spectrum. On the whole, however, 
systematic differences in the absolute magnitudes arising from personal judgments 
of the two observers appear to be certainly less than o™1. 

4.3. However, a wider investigation is needed. Not only have there been 
two different observers; there have been two telescopes, two spectrographs, 
different photographic emulsions differently calibrated for photometry, two micro- 
photometers, etc. Originally there were only 5 stars on the Pretoria list (Table III: 
Nos. 94, 95, 98, 100 and 106) for which absolute magnitudes had already been 
obtained at Victoria, and these covered a range of 2-5 magnitudes only, while the 
Victoria calibration-curve covers altogether 6-7 magnitudes. Further examina- 
tion of material available at Pretoria and Victoria revealed a much more satisfactory 
overlap. There were 39 stars outside the original programme for which radial 
velocity work had been done on photometrically calibrated plates both at Victoria 
and Pretoria; to help widen the range of intercomparison, a number of additional 
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spectrograms were specially obtained at Pretoria. Eventually there were thus 
44 stars in common, with 128 spectrograms from Pretoria and 88 from Victoria. 

Thanks to the kindness of Dr Petrie, all the microphotometer work on these 
plates was carried out at the Dominion Astrophysical Observatory, following the 
procedure he used in his previous work on northern stars. His comparison of the 


TABLE I 
Comparison of the equivalent widths of Hy between Pretoria (R) and Victoria (DAC) plates 


R DAC 

Star W n W n 

I HD 23408 6:10 3 6°70 4 
2 23793 6:90 4 7°20 2 
3 30677 2°12 2 2°93 3 
+ 35730 5°92 I TIS 2 
5 36013 6:18 3 5°98 2 
6 36895 6:06 2 6°32 I 
4 36960 3°28 2 4°14 2 
8 37018 4°13 2 4°01 I 
9 37128 1°98 6 1°42 3 
10 37209 4°53 6 4°95 4 
II 37742 1°50 I 1°74 2 
12 46056 2°03 2 3°80 2 
13 46202 3°82 2 3°96 2 
14 46573 2°17 3 2°36 2 
15 46883 2°33 I 3°34 2 
16 50868" 4°01 2 5°54 2 
17 66665 3°17 3 3°36 2 
18 138485 5°54 9 5°90 I 
19 142096 7°13 5 7°87 I 
20 144218 ~ 6:29 a 5°64 2 
21 149438 3°62 8 3°85 4 
22 154445 4°58 2 4°04 2 
23 160233 3°89 3 4°62 2 
24 162365 4°81 2 5°63 I 
25 163472 4°80 4 4°34 2 
26 167785 5°98 2 6°10 2 
27 170580 641 3 5°44 2 
28 174391 5°13 4 5°62 2 
29 175803 5°98 2 6°80 2 
30 176304 4°42 3 4°39 2 
31 176542 6°42 2 7°86 I 
32 178129 2:29 3 2°13 2 
33 181858 6°41 I 6°58 2 
34 181963 4°76 I 4°82 2 
35 184279 3°41 3 3°70 2 
36 187320 3°52 I 2°98 2 
37 189550 4°54 2 5°30 2 
38 189921 71S 3 6°62 2 
39 196025 5°16 4 5°29 2 
40 201254 5°86 2 5°52 I 
41 207563 5°56 4 5°67 3 
42 209008 5°69 3 6°35 2 
43 215733 2°64 2 1°99 2 
44 219188 2°84 2 2°69 2 
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Fic. 3.—Comparisoz of Hy intensities from spectra obtained at Pretoria (Radcliffe Observatory) 
and Victoria (Dominion Astrophysical Observatory). Open circles have a smaller weight (stars 
with one plate only). 


individual Hy intensities is plotted in Fig. 3, with Victoria values as abscissae and 

Pretoria as ordinates. ‘Table I shows the results for Hy, and Table II for A4388A, 
Her. ‘The computation of regression lines by least squares gave 

Wpro= —0°03 + 1-044. Wp (5) 

Wp= +0°06 + 0°938. Wp ao (6) 

Taking these relations at their face value, the Pretoria equivalent widths are about 


5 per cent smaller than the Victoria equivalent widths. In view of the general level 
of errors involved, this difference has been neglected in what follows. The 
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Victoria calibration curve and the Victoria spectral-type corrections have been used 
0) _ unchanged for the Pretoria measures, and it should be possible to regard the 
, Pretoria values as on exactly the same footing as the Victoria values. 

5. Absolute magnitudes of southern B-stars.—The following data are tabulated 
in Table III which comprises all stars investigated in the present Hy-programme. 


Column : Column : 
1 Current number g Colour index, B—V. 
2 HD number 10 3x colour-excess, 3CE. 
3 Galactic longitude, /! 11 Number of spectra measured, 1. 
4 Galactic latitude, b!. 12 Equivalent width, W. 
5 Spectral and luminosity classification, 13 Standard error of W in units of 
Sp. oor A, + AW. 
6 Proper motion in right ascension,u,. 14 Absolute magnitude, M. 
7 Proper motion in declination, 1 ;. 15 Standard error of Min units of 
8 Visual magnitude, V. o™-o1, +AM. 


16 Distance of star in kiloparsecs, r. 


6. The luminosity-spectrum array.—Fig. 4 (a) and (b) represents the luminosity- 
spectrum relation for all the 137 stars listed in Table III, which may be compared 
(in Fig. 4 (b)) with the corresponding relationship for 313 northern stars (derived 
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from Petrie, 1956b, p. 299; and 1958a). In obtaining mean magnitudes for 
each spectral sub-type of our southern stars, (Tables IV a and IVb), we have 
followed Petrie’s procedure of excluding, more or less arbitrarily, certain thinly 
scattered high-luminosity stars, regarded as supergiants. Fig. 4(a) shows the 
individual southern stars, while in Fig. 4 (b) these have been brought together into 
six groups, with the addition of the corresponding mean values for the northern 
stars. It is apparent that southern and northern stars show very much the same 
relation between mean absolute magnitude and spectral type. It should be noted 
that on the whole the northern stars are nearer the Sun than the southern stars, 
but that the distribution in absolute magnitude appears to be essentially the same 
for the greater volume and depth of space now sampled. 


TaB_e [Va 
Sp. M n 
O 6 — 4°95 4 
O 65 ~$s I 
O 75 — 5°35 4 
O 8 — 4°35 2 
O 85 — 5°80 2 
O 9 — 4°91 8 
O 9°5 — 4°48 14 
Bo —4°17 II 
B os — 4°30 15 
B 1 — 3°16 20 
B is —2°2 I 
B 2 — 1°89 25 
B 3 — 1°62 13 
B 5 —I°1o 2 
B 7 —1'5 I 
TaB_eE IVb 

Mean Mean Spectral 

type M range n 
O 72 =5°n5 06-0 8&5 13 
O 93 — 4°64 O09 +095 22 
B o3 — 4°25 Bo+B o5 26 
B ro — 3°02 Bis+8 1 21 
B 2:0 — 1°89 B 2 25 
B 3°5 — 1°46 B3-87 16 


7. A uniform system of southern star magnitudes and colour indices. 

7.1. To obtain the distances of the stars with which we are dealing, we need to 
know both the apparent visual magnitude V and the interstellar absorption. The 
absorption has to be estimated from the colour excess, so that we need to know the 
star colours. Examination of the available data revealed a rather heterogeneous 
material of magnitudes and colours on different ‘‘systems’’. An attempt has 
been made to reduce all of these to acommon svstem. The majority of the stars 
have Cape magnitudes, and reduction formulae to adjust Cape measures to the 
(B—V)-system have been given by Stoy, so that what may be called the ‘‘ Cape 
(B—V)-system’”’ has been used as a standard in the present paper. The colour 
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excess has been taken to be the difference between the (B— V)-value, found in this 
way, and the unreddened (B—V), given by Johnson (1958) for spectral types 
O8 to B7 and three luminosity groups (lab, II, and III-V). The interstellar 
absorption for the V-magnitudes has been assumed to be 3-0 x CE, i.e., three times 
the colour excess, following discussions by Johnson—Morgan (1953), Morgan- 
Harris—Johnson (1953) and Morgan (1958). 

Entering the calibration-formula with the measured W, we obtain the absolute 
magnitude M’; the addition of the spectral-type correction AM’ gave M in 
Column 14 of Table III, which together with Column 8 and Column 10 led to the 
distances r (in kiloparsecs) of Column 16: 


r=10~* [antilog o-2 (V— M’—AM’+5-3CE)]. (7) 


(1) Colour indices 
7.2. The available sources of colour indices for our stars are: Stebbins 
(1940); Oosterhoff (1951); Schilt (1952); Sharpless (1952); Morgan, Code 
and Whitford (1955); Johnson (1955); Hiltner (1956); Stoy (1958a, 1960a); 
Willstrop (1959). 
7.2.1: Stoy.—The Cape colour-indices SCI may be converted to (B—V) 
by a relation given by Stoy (1958 a): 
B-—V=0°85 SCI + 0:26. (8) 
7.2.2: Morgan—Code—Whitford.— A list by Stoy (1958b) has 12 stars in 
common with Morgan—Code-Whitford (1955); Fig. 5 shows the MCW colour- 
index C, as a function of Stoy’s (B—V). The relation is very smooth and can be 





C,(McwW) 
! 











wh 
m m m 
0-0 +05 +10 


(B-V)C 
Fic. 5.—Relation between Morgan—Code—Whitford’s C, (MCW) and the Cape values (B-V)C. 


used to convert all the relevant MCW-colours to (B—V). Furthermore, this 
curve agrees closely with a formula given by Morgan—Harris—Johnson (1953, 
equations (2) and (3)) for the relation between C, and (B—V); they found: 

C, = —0°144+0°483 (B—V). 
From this we may write 


B-V=2:07C,+0°30. (9) 
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The mean of the (B — V)-values resulting from Fig. 5 and from the above expression 
was used to convert C, to (B-—V). Data were obtained for 29 stars from this 
source. 

7.2.3: Stebbins.—The Stebbins-values C, have been regarded as the same as 
the MCW-values C,. 

7.2.4: Oosterhoff.—Colour indices measured by Oosterhoff, which he calls C, 
and C,, may be converted to (B— V) toa first approximation by 


(B-V) = ~32 [(0-571 +0565 C1) + (0367 +0797C,)] +016 (10) 


This formula results from a combination of the following set of relationships 
involved 


Cy’ =0°571 +0°565 C, (11) 
C,’ =0-367+0-797C, (12) Oosterhoff, 1951 
C*=0'5 (Cs’ + Cs”) (13) 
ICI =1-50C* (14) Stebbins, 1940 
ICI =0-93 SCI +0-11 (15) 

B—V=08s5 SCI +0-26 (16){ Sty 79580. 


Small corrections to the above formula (10) for (B—V) were evaluated by 
using 87 Oosterhoff-stars for a comparison with (B— V)-values deduced from the 
other sources mentioned. In this way adjusted colour-indices were, in addition, 
obtained for those 22 of the above Oosterhoff-stars for which there are no other 
observational values. 

Stoy (1960b) has reduced Oosterhoff’s indices C, and C, by the expressions 


(B—V),=1°042C,+1°313, and (17) 
(B—V),= 1-481 C, + 0°948, (18) 
which are based on the comparison of Oosterhoff’s colours with the colours for 
18 similar reddened B-type stars. ‘The mean of these (B—V),-and (B—V),- 


values makes those of Oosterhoff’s stars whose colour-indices exceed +0™-05 
by about o™-04 bluer than they are with our above reduction. 


(II) V-magnitudes. 

7.3. After a good deal of intercomparison among themselves and with Cape 
(based on 80 southern stars) it was decided that the expressions listed below 
(equations (19), (20), (21), (22)) give the best reduction to a self-consistent system 
(Cape- B—V) that can be obtained from data available at present. This step 
required certain additional empirical corrections to the various conversion- 
formulae; the latter can be deduced from the references given below (Sections 


7.3.1-4). 
7.3.1: Stoy (1958a): 
V (Cape) = 1Pv = SPv + 0-08 SCI — 0:06. (19) 
7.3.2: Morgan-Code—Whitford (1955, p. 42): 
V (Cape) = [magnitude in MCW, l.c., Table II] — 0-12. (20) 
7.3.3: Stebbins (1940, p. 24): 
V (Cape) = m (Stebbins) + 0-05. (21) 


7.3.4: Oosterhoff (1951, p. 320): 


V (Cape) =0-951 m (Oosterhoff) — 0-28 C,+ 0°15. (22) 
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7.3.5: Magnitudes given by the other authors mentioned in Section 7:1 
are on the (B—V)-system and need no correction. Where a star has been 
measured by more than one investigator, a straight mean has been taken of the 
available V-values, if necessary corrected as just described (‘Table III, Column 8). 
In the same way, straight means have also been used to combine (in Column 9g) 
the adjusted colour-indices derived from different sources. 


8. Calibration tests with mean proper-motion parallaxes and galactic rotation. 
We have seen that, where the measures overlap, the Pretoria and Victoria results 
agree quite closely. Petrie (1953, 1958b, 1960) has devoted much care to obtain- 
ing as reliable a calibration of the (Hy, /)-relation as is possible with our present 
knowledge of the stars accessible to northern observers. It does not appear 
practicable yet to add much calibration evidence from the southern sky, but a 
number of checks have been made and the results are summarised here. 

8.1. Clusters.—Stars Nos. 94, 95, 98, 100 and 106 are fairly certain members of 
the Scorpius-Centaurus cluster. Blaauw’s (1946, see also 1958) cluster-parallaxes 
together with the V- and (B— V)-values, lead to a mean absolute magnitude 


Mo=—3™0. Themeanfrom Hy is Mg= —2™1. 


8.2. Individual parallaxes.—Table III includes 14 stars for which there exist 
other parallaxes (trigonometric, dynamic, etc.). ‘The mean of the differences of 
the corresponding Hy-parallaxes minus the other values is Az = — 0”-o01 + 0”-001. 
The similar result for 53 northern stars was + 0”-002 + 0”-004 (Petrie 1953, p. 265). 

8.3. Proper motions.—We have followed the procedure that Petrie (1956b) 
used for the northern stars. Proper motions for the southern stars were taken 
from the Cape Annals, the Yale Transactions and the General Catalogue. Of 132 
stars with known colour excesses, there were gg stars for which p, and p, were also 
available. Although not used here, new Cape proper-motion data for 11 mors 
became available later (Stoy 1960 a) and are included in Table III. 

8.3.1: The v-components of the proper motions were obtained using the 
convenient diagrams devised by Pearce (1927), which assume the solar apex to be 
at «=271°, d5= +28° (1900). Then 

- _k>vsina 
70 Vy Ssind’ (23) 
where k=4°738, 
V,,=20°0 km/sec, (solar velocity), and 
A= distance from star to apex. 


The stars are arbitrarily separated into three groups (although this makes each 
group undesirably small), viz., 


*“ dwarfs ’”’ with Msgp> —2™-0 
** intermediates ”’ with —2™-0> Msp> — 4™0 
“giants ”’ with Msp < —4™-0 


The frequency-distributions of our 19 dwarfs, 32 intermediates, and 48 giants are 
shown in Fig.6{a). ‘The comparison with the northern stars in Fig. 6 (6) demon- 
strates the very different distribution of the northern material (Petrie 1956a, 
1956b, 1958): the giants in the latter are concentrated at about 1000 parsecs, 
while for the southern stars the maximum lies at about 2700 parsecs. 
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. 8.3.2: We need to estimate the standard errors of both the means of the Hy 
parallax 7sp, and of the mean proper motion parallax 7,. 
ach Since 
) M=m+5+5logiy7, 
we have 
AM = 5 logy) e (Am/z) = 2°17 (Az) /z. (24) 
| Petrie (1953, p. 265), on the basis of internal agreement of his measures, gives ‘‘ as 
sare a generous estimate of the probable error of a single spectroscopic absolute 
1on- | Magnitude’? +o™-33; or,say, +o™-5 asthe standarderror. The corresponding 
56a, ff parallax-error would be, therefore, 
secs, ff O° 
: q Asp = = = Tsp >= a O°2 Tsp» 
2°2 
15* 
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or, of the mean parallax of n stars: 
o°-2 Tsp 
+——. 2 
* (25) 
This takes no account of errors in the correction for interstellar absorption, so that 
the total Azsp are probably somewhat greater than estimated here. 
For the error of the mean proper-motion parallax we accept the validity of 


where R is the mean linear random motion in one coordinate, taken to be 5 km/sec. 
The mean 1/sinA is assumed to be 1-75, and for ‘‘intermediate’’ stars the mean 
Ap/v is 1-2, as found from the General Catalogue. The mean value of the quantity 
under the root is then 1-3. 

8.3.3: We thus obtain: 


A7sp = 





Az, (intermediates) = + 2:07,/4/n, (27) 
Az, (dwarfs) = +1-07,/+/n, (28) 
Az, (giants) = +3:87,/+/n. (29) 


We then find for 19 ‘‘ dwarfs ’’: 
7, =0":0033 + 0”-0009 
Tsp = 0":0039 + 0”-0002, 
for 32 ‘‘ intermediates ’’: 
@,,=0"-0005 + 0”-0002 
Tsp =0"-O0015 + 0”-0001, 
and for 48 “‘ giants’: 
7, =0"-0008 + 0”-0004 
Tsp =0":0005 + 0”-0000. 

8.4. Galactic rotation.—The spectroscopic distances based on Hy-intensities 
may also be compared with distances deduced from the galactic-rotation effect, 
again following the procedure used by Petrie for northern stars. Petrie calculated 
a mean distance r, from the formula 
se 1 Dpsin2L cos*b 

9 A sin? 2L cos*d ’ 
where p is the observed radial velocity of each star after correction for the solar 
motion and the K-term; A is Oort’s constant, taken to be 17-5 km/sec/kpe, in 
accordance with a’ recent discussion of the B-star material by Feast and 
Thackeray (1958; see also Petrie 1956b); L is the galactic longitude of the stars 
measured from the longitude of the centre of the Galaxy; b is the galactic latitude ; 
the distance r is expressed in kiloparsecs. 

8.4.1.—The standard error of the spectroscopic distance, by analogy with 
equation (25), may be taken to be 





se o'27 
Afsp = + ~~ y (30) 


The standard error of the mean galactic-rotation distance may then be 
estimated as 





Ar, += s (parsecs), (31) 


which follows from the assumption that the random velocity of the average star 
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plus its observational error add up toa radial velocity AV =6 km/sec; furthermore 
that p has on the average half of that galactic-rotation amplitude which would 
correspond to the distance r, that is, p= 0-00875r parsecs. We can then calculate 
the above expression (31) from 
Ar, = + at ve. (32) 
8.4.2.—There is, however, a difficulty. ‘The southern stars are on the whole 
more distant than the northern; thus, the approximate galactic-rotation formula, 
which assumes that r is small relative to the Sun’s distance from the galactic 
centre, breaks down for roughly half of them. ‘This is unfortunate, for it is 
precisely with the ‘‘ giants ’’ which constitute the more distant half of our material 
that we most need to test the Hy-results. 

Confining ourselves to the approximate formula, we consider only stars not 
more than 1250 parsecs from us. We divide these into two groups, according to 
their Hy-distances : 

(A) ... 250<rsp < 750 parsecs (25 stars): 
7g=459 + 144 pc 
7sp= 386+ 16pc 
(B) ... 750<rsp < 1250 parsecs (25 stars): 
7g=1561 +144 
Tsp=1045+ 45. 

Both tests, that with proper-motions (Section 8.3.3) and that with galactic- 
rotation distances, refer to too few stars for the results to have a decisive weight. 
There is a tendency for the mean Hy-parallaxesto be larger (or the mean 
Hy-distances smaller) than those given by proper motions or galactic rotation. 
The difference is not significant for the nearer stars, mainly those we have called 
‘“‘dwarfs’’. It may be just significant for the ‘‘intermediates’’, while about 
the ‘‘giants’’ we can say nothing. 

8.4.3.—If we now attempt to use instead of the velocity-displacement of the 
usual stellar lines, that of the interstellar K-line, and then deduce ‘‘ calcium 
distances ’’ by means of the galactic-rotation formula, we find 





Group A (13 stars) Group B (22 stars) 
calcium 7g = 525 + 190 pe 1163 + 148 pe 
stars fsp=425+ 27pec 1044+ 45pc 


We note that the r-values agree reasonably well within each of these groups. 
This finding, however, is completely in conflict with the expectation that ‘‘ effec- 
tively ’’ the calcium should be only at half the distance of the star in whose spec- 
trum it is observed. As far as the northern stars are concerned, Petrie’s (1956b, 
p. 295) data seemed to confirm this expectation. Although the number of 
southern stars available for this test is certainly small, this of itself seems an 
inadequate explanation. The above figures clearly indicate that here, in our 
southern material, the calcium appears to be located at about the same distance as 
the stars. Here we note that the considerable number of stars involved in an 
Oort-type analysis by Feast and Thackeray (1958, p. 133) did not exhibit such a 
discrepancy. 

There are two remarks which may be made. First, high-resolution measures 
of the interstellar calcium lines by Adams and others showed that these lines are 
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almost always multiple, the calcium being concentrated into discrete clouds along 
the line of sight. Second, while customarily we think of these clouds as being 
fairly uniformly distributed near the plane of the Milky Way, so that the line of 
sight intersects a number of clouds of the order of ten per kiloparsec, it seems highly 
likely that there is some systematic distribution of cloud-frequency related to the 
spiral arms of the Milky Way. The view is held by some that the Sun is near the 
inner boundary of a spiral arm. In that case the first kiloparsec or so along the 
line of sight in the Milky Way north of the equator will be more or less within a 
spiral arm, while in the southern Milky Way it will mostly run through an inter- 
arm region. If calcium clouds are a part of the arm population, but are compara- 
tively rarely found between the arms, the results just quoted become a little less 
puzzling. But in any case, although selection effects cannot be ruled out, we 
feel inclined to reject the calcitum-method as a safe guide to the distances of the 
southern B-stars. 

8.4.4: There is one further point about the B-star distances that must not be 
forgotten, and that is the uncertainty of the correction for interstellar absorption: 
there are the ordinary errors in (B— V) and, in addition, the errors in the assump- 
tion that the absorption corrections for V-magnitudes is 3-0-times the colour 
excess in (B—V). In conclusion, we note that the average correction 3CE in 
Table III is o™-71 for ‘‘dwarfs’’, 1™-o1 for ‘‘intermediates’’, and 1™-40 for 
“giants ’’. 
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Fic. 7.—The colour-excess array (both hemispheres). 
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g. Relation of results to the structure of the Milky Way. 

g.1.—It is of considerable interest to try to fit the results of the preceding 
sections into our present picture of the galactic structure in regions within a few 
kiloparsecs of the Sun. For this part of the discussion we have also included 
material from the northern hemisphere, namely 405 O- and B-stars for which 
Petrie (1953, 1956a, 1956b, 1958) has obtained Hy-luminosities. To derive 
distances for these, all published magnitudes and colour indices were re-reduced 
to the (B—V)-system as uniformly as possible, whereby the derivation of the 
corrections followed a procedure similar to that used for the south. In particular, 
comparisons of Stebbins- and MCW-data with Hiltner’s (1956) list of (B—V)- 
measurements could be based on 48 stars for V-magnitudes, and on 51 stars for 
(B—V)-colours. The details will be omitted here. 


145 





Fic. 8.—Colour-excess diagram of all B-stars with measured Hy-luminosities in both hemispheres, 
out to nearly 4 kiloparsecs. The diagram represents the distance at which the star is affected by a 
certain absorption, so that where points are far from the centre, the general absorption is small, and 
vice versa. The three broken circles correspond to distances from the Sun of 1, 2 and 3 kiloparsecs, 
respectively. The outermost closed “‘ curve ’’ represents the relation between colour-excess and dis- 
tance for CE=o™-8, In addition short extreme “ arcs”? for o™-g, 1™-o and 1™-1 are indicated in 
the 60° galactic longitude sector between 345° and 45°. All the iso-absorption lines are drawn at 
CE-intervals of o™-1. 


g.2.—If colour excesses, for both northern and southern stars, are plotted 
against distance (Fig. 7), a steep increase of reddening appears up to a distance of 
about 1 kiloparsec, but comparatively little more thereafter. ‘The chief cause of 
this is probably observational selection. For example, Bo stars beyond 1 kpc 
will mostly appear fainter than 9™ and may have the apparent colour of a G- 
or K-type star so that many, or possibly most, of those within 3 kpc have not yet 
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appeared in radial-velocity or photometric programmes. ‘There certainly has 
been selection in favour of apparently blue stars lying in relatively transparent 
regions, and against those subject to unusually heavy reddening and general 
absorption. To some extent even stars within 1 kpc must have been affected by 
this selection. 

9.3.—Taking our results at their face value, we have plotted a series of contours 
of equal colour-excesses (in steps of o™-1) in the plane of the Milky Way (Fig. 8). 
The interstellar absorption affecting our OB-stars appears to be greater in the 
region 330° </! < 165°, which lies mostly in the northern sky, than in the remainder 
of the Milky Way, 165° </'< 330°, lying mostly south of the equator. Not too 
much weight should be attached to this until the question of differences in observa- 
tional selection between north and south has been examined more thoroughly. 
Here we should also note that, for instance, our data do not even cover all the 
already known O-B3 stars brighter than 7™-5 in the southern sky. 


228 218 208 158 
—— ' 























Zz 
CS 278 








Fic. 9.—The B-stars in the sector 1} =275°+: 
©O=stars with 3CE <1™-5. 
© =stars with 3CE>1™+s. 
Lightly dotted shading= H1 areas with densities d=0-6 atoms|cc. 
Heavily dotted shading = H1 areas with densities d=1 atom|cc. 


The large unshaded areas have a density of only 0:2 atoms/cc,—except for an extended region 
enclosed by the curve on the top left, the smaller elliptical area underneath, and the three very small 
enclosed regions; for all these five last-mentioned areas the density is only 0-05 atoms|cc. 


9.4.—In Figs. 9, 10 and 11 we have plotted, with the Sun as centre, B-star 
distances as derived from Hy measurements. Fig. 9 represents the B-stars in the 
particularly interesting sector around /!= 275° (see Section 9-7); Fig. 10 reaches 
to r=1'5 kiloparsecs, while Fig. 11 combines all the 465 measured B-stars with 
known 3CE up to 7kiloparsecs. A number of characteristic constituents of the 
Milky Way are added, as described further below in Section g.5 and in the captions 
to the figures. It should be noted that of the 333 northern stars with known 3CE 
only 17 haver>1-5kpc. Furthermore, of the 316 northern stars within 1-5 kpc 
there are 22 stars with 1-5<3CE<2:o and 12 stars with 3CE>2-0. This 
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Fic. 10.—Galactic objects within 1:5 kiloparsecs. 


(B-stars, H t-areas, H t1-regions, clusters, associations, spiral-arm patterns.) 
@ =B-stars, both hemispheres, with 3CE <1°5. 
© =B-stars, both hemispheres, with 1-5 <3CE < 2:0. 
© = B-stars, both hemispheres, with 3CE > 2:0. 
/\=clusters. 
A=associations. 
Shaded portions=H t-areas with densities d>1 atom|cc. (The adjacent dotted lines 
enclose areas of density 0-6 and less.) 

©O=H ti-regions. 

‘--'-- *——-+-—-+=Becker’s spiral-arms. 

Diagonal hatching= Elsdsser-Haug’s spiral-arms. 
----------- = Mills’ spiral. 
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Fic. 11.—Galactic objects out to 7 kiloparsecs from the Sun; (465 B-stars, 82 clusters, 55 OB- 
associations, 36 H 11-regions): ( 
@ =southern B-stars with 3 CE <1°5. 
© =southern B-stars with 1-5 <3CE <2:0. t 
© =southern B-stars with 3CE>2:0. I 
@ =northern B-stars with r<1-5 kiloparsecs and with 3CE <2:0. ] 
%*=northern B-stars with r<1-5 kiloparsecs and with 3CE>2°0. 
4 =northern B-stars with r>1°5 kiloparsecs and with 3CE <1°5. 
Xx =northern B-stars with r>1°5 kiloparsecs and with 1-5 <3CE < 2:0. ] 
%*=northern B-stars with r>1°5 kiloparsecs and with 3CE>2°0. 
Shaded portions=H t-arms (density d> 0-6 atoms|cc). 
/A\=clusters. § 
A=associations (j—> =Dieter-distances). ¢ 
© =H u-regions. 
tem em meee = Becker-arms. t 
Diagonal hatching = Elsdsser-Haug-arms. 
---------- = Mills’ spiral. t 
Note—The star No. 40 (l=257, r=5-°78) should be added to Fig. 11 above. 
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distance-range (i.e., Fig. 10) contains, in addition, 44 clusters, 33 OB-associations, 
and 21 H 11 regions from our lists (Section 9.5). 

When we come to trying to locate our B-stars in supposed spiral arms, we 
find ourselves faced by a complicated problem and insufficiently supplied with 
reliable information about the spiral arms, both with respect to the radio- as well as 
to the optical evidence. 

In the first place we have the well-known maps of the neutral hydrogen distri- 
bution in the galactic plane, based on 21cm measures. There is a certain amount 
of ambiguity in the final interpretation of the observations, and in that part of the 
sky in which we are chiefly interested here, south of the equator, Australian 
workers are just revising and improving their results, introducing a different 
galactic model which involves a specified expansion-motion of the hydrogen gas 
(Kerr 1961). I am very grateful to Dr F. J. Kerr for this information and for his 
permission to use some details of his new contour-map in the present discussion 
of our new B-stars positions, the result of which he considers as ‘‘a first clue’’ 
in view of the fact that there has been little success in previous attempts to relate 
optical- and radio-spiral-features. 

On the one hand, it is certainly questionable whether spatial structure marked 
by aconcentration of neutral hydrogen is to be identified with spiral arms consisting 
of a concentration of O- and B-stars, gaseous nebulae, H II regions, etc. 

On the other hand, a considerable body of experimental, i.e. observational, 
evidence has accumulated concerning a possible spatial connection between the 
H1 ‘‘arms’’, the galactic-cluster arms according to W. Becker, the Elsasser-Haug 
spiral-structure, the Mills-spiral, and the distribution of the southern and 
northern B-stars, open cluster, OB-associations and H 11 regions. It is profitable 
at this stage to attempt a survey of the available evidence. 

g.5.—Fig. 11 brings together in one diagram all the above-mentioned objects. 
Taking this diagram at its face value, it becomes, in particular, at once evident 
that the stellar objects and the Hi regions appear to be situated between the 
H I-arms. 

Fig. 11 assembles: 

(a) Southern B stars (this paper, Table III). 

(6) Northern B stars, deduced from Hy-measurements given by Petrie 
(1953, 1956a, 1956b, 1958). 

(c) Combined and revised H 1-contours from both hemispheres, according to 
the 21cm work of the Leiden and Sydney observers; contours restricted to 
regions denser than, or equal to, o-6atoms/cc; after Kerr, 1961; see also Oort- 
Kerr-Westerhout, 1957, 1958a, b; van de Hulst, 1958; Bok, 1959; Struve, 1959). 

(d) Hit-regions (36); after Morgan-Sharpless-Osterbrock (1952); see 
Becker (1961); also Bok (1959). 

(e) galactic clusters (82); after W. Becker (1961; see also 1956); three 
spiral arms, 1-5 kiloparsecs apart, are suggested by these clusters over an arm-length 
of 4-5 kpc each (centrearm: Orion; south: Carina, Scutum; north: Perseus) ; 
thickness of Orion arm about 1-5 kpc near Sun. 

(f) spiral arms indicated by Elsisser-Haug’s (1960) surface-photometry of 
the Milky Way; thickness of each arm about 0-6 kpc. 

(g) spiral arms proposed by Mills (1959) from his observations of the radio- 
continuum radiation; this structure was suggested as mainly to represent the 
distribution of ionized hydrogen. 
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(h) OB-associations, after Blaauw, Code, Johnson, Greenchik, Gulledge, 
Hardie, Hiltner, Morgan, Schulte, Seyfert, Whitford. For references see 
Becker (1961) and Dieter (1960); these two sources, together with additional 
data compiled by van de Hulst (1958) and Bok (1959), supplied 55 associations. 

(i) OB-association-distances discussed by Dieter (1960); positions found 
photometrically and those derived from galactic-rotation are joined in Fig. 11 
by straight lines with arrowheads. 

g.6.—Examining the constituents of this diagram in turn, we note that the 
groups (a) and (b) clearly fall outside (c) and (g), but appear well-connected with 
the groups (4), (e), (f), (A), (i): 

g.7.—A conspicuous branch of B-stars, especially reddened ones, stretches in 
the galactic direction /'=275° +, which lies well between two prominent HI arms 
(Fig. 9). ‘These stars appear to be of particular significance : 

If we examine the longitude-colour distribution of all the 132 stars with known 
colours, and denote with r the reddened stars with 1-5 < 3CE < 2:0, and with wr the 
very reddened stars with 3CE > 2:0, we obtain the following distribution, referred 
to five different ranges in galactic longitude, the last line giving the percentage 
of the (r + vr)-stars in the sector : 





o —60 130-190 220-270 270+280 280-360 
17 stars 15 38 30 32 
Ir 2r 6r + 377 or + 6cr 4r+5uvr 
6 13 24 50 28 per cent 


(None of the stars in our Table III falls into the longitude-ranges 60°—130° and 190°-220‘). 


The conspicuous concentration of the r- and vr-stars between 270° and 280° 
becomes even more emphasized if we widen this range somewhat to cover the 
sector 265 °-287°, centred around 275°. Here we have now 34 stars, of which no 
less than 10 are r-, and 9g vr-stars, i.e., together 56 per cent of all the stars in this 
sector. 

We compare this with the material in all the other sectors: there remain 98 
stars, comprising only 12,7 as well as 5 vr-stars, i.e., together only 17 per cent of the 
total. ‘Therefore, we can state that a sector of less than twenty degrees around 
l'=275°—which lies essentially between the ‘‘hydrogen arms’’—contains 
(considered in percentages) no less than 43 times more ‘‘ very reddened ’’ stars, 
and 24 times more of the ordinary ‘‘ reddened ”’ stars, than there are in the whole of 
all the other sections of our southern material. 

We may also note other important regions outside the H I-areas, situated in 
northern sectors, particularly those centred about /'= 40° + and/'=go°+. 

The 275° +-region is shown in detail in Fig. 9, where the newly-determined 
B-star positions are compared with neutral hydrogen contours. Most of the stars 
are seen to appear in a region well outside the Sagittarius arm. Kerr’s successful 
recent Milky Way model, involving, as it does, a combination of two motions,—an 
outward-flow of the gas from, and its rotation about, the galactic centre—is 
clearly in general qualitative agreement with the distribution of stars shown in this 
diagram. Kerr has characterized the conclusions from Fig. g by suggesting that 
if the material of the arms is moving outwards from the central regions of the Milky 
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Way, the outer part of each arm should contain the oldest gas, and one may therefore 
conclude that this part of the gas may have been largely condensed by now, i.e., 
used up to produce stars in the outer part of the arm,— indeed just as indicated by 
the data shown in Fig. 9. Moreover, it is apparent that any given general H1 
pattern should, under the influence of these combined motions, have a relatively 
short life-time, and stars and hydrogen should dissociate from each other at early 
evolutionary stages. 

g.8.—In connection with group (d) the extensive recent Mount Stromlo work 
on the southern H1ir regions (Rodgers, Campbell and Whiteoak, 1960) is of 
special importance. Dr Rodgers kindly showed to the author a relevant sheet of 
the Canberra H«-atlas now in preparation: it demonstrated that our ‘‘ reddened 
sector’ around 275° is particularly rich in markedly small H 11-regions; more of 
these appear here per unit area than anywhere else. This points to the radiation 
of numerous OB-stars and, because of the pronounced reddening-effects observed 
in this sector, to the presence of strong absorption by dense interstellar dust. 

9.9.—As to (f ), this recent optical work by Elsisser and Haug (1960) is based 
on a two-colour photoelectric photometry of the Milky Way. From these mea- 
surements the authors have deduced outlines of three spiral arms which show little 
detailed agreement with the H1 arms given by the 21cm work. Plotting histo- 
grams of frequencies for our B-star material, however, we find that these stars show 
some tendency to lie within the Elsisser-Haug arms. We realize, of course, that 
these spiral arms should at present merely be considered as a schematized picture. 

g.10.—Of particular significance for our problem is W. Becker’s (1961) 
extensive analysis of 82 open clusters. Here, it is pointed out that there is no coin- 
cidence between the positions of these clusters and the radio-astronomically 
determined H1 regions but, on the contrary, a fundamental difference in their 
distributions. Becker’s method of the photometric determination of cluster- 
distances is estimated to be correct within 10 per cent of the distance, and there 
seems no possibility of reconciling the cluster positions and the 21 cm results. As 
Becker puts it: ‘‘the question arises whether a plausible change in the law (of 
galactic rotation) might lead to a coincidence of the distributions ’’. 

g.11.—A similar result can be derived from a recent analysis of the OB-associa- 
tions: in comparing distances of associations obtained, first by photometry, and 
secondly from galactic-rotation theory using the observed radial velocities, Dieter 
(1960) found in most cases appreciable discrepancies. ‘The conclusion was drawn 
that the ‘‘ application of such a simple (circular-motion) model to the determina- 
tion of the distances of neutral hydrogen clouds without independent distance- 
measurements may lead to serious errors ’’. 

It would appear, therefore, that H 1-positions in our Fig. 11 may need to be 
modified. We recognize a suggestion of such a correction indicated in Dieter’s 
Fig. 3 (1960, p. 55), the relevant feature of which is added to our Fig. 11; here, 
the heads of the arrows denote the galactic-rotation positions, and their tails the 
corresponding photometric positions. If we assume the photometric distances 
to be correct, we may apply these vectors (reversed) as tentative corrections to the 
H1 positions: if we do so, the latter ones come to a large extent in coincidence with 
the plot of stellar objects. Manifestly, much future work is needed. 

Finally, it is interesting to note that Dieter found evidence for a coincidence 
with strong H1-concentrations, situated within the accepted 21 cm arms, for 31 
out of the 40 investigated OB-associations. 
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g.12.—The whole subject of galactic rotation, hydrogen distribution, lumino- 
sities and distances appears today more fluid than ever. As a postscript we may 
mention here another recent paper: Thiessen (1961), investigating the origin of 
the polarization of star-light, comes to the conclusion that stellar radiation appears 
to be composed of unpolarized thermal radiation and polarized synchrotron 
radiation, and that the observed phenomenon is a stellar rather than an interstellar 
one. As a consequence, Thiessen points to possible implications of the 
proposed new mechanism, which may require that, e.g., absolute magnitudes, 
intrinsic colours, colour-magnitude diagrams, etc., may ‘“‘ partially have to be 
corrected’. We also mention important discussions by Johnson and Iriarte 
(1958). 

10. Possibilities of photoelectric work.—Recent investigations carried out at the 
f/18 coudé focus of the 36-inch reflector at the Cambridge Observatories, in which 
the integrated intensities of selected features of stellar spectra are measured 
directly, have clearly shown great potentialities. Following the construction of a 
new spectrometer, R. F. Griffin (1960) made a few exploratory measures during 
one night in June 1960 of Hy in B-type stars, which showed that it would be 
possible, with the present instrumental arrangement, to carry out a programme 
down to 7™5. 

Subsequently some 1200 observations were made by the writer on 18 nights 
during 1960 October and November, giving Hy-intensities for 187 B-stars. ‘The 
measurements were made with two photomultipliers, one receiving light from the 
spectrum range AA 4325-4357 A, the other from two comparison regions, AA 4255- 
4289A and 4393-4427A. The general procedure was that followed in earlier 
work by Griffin and Redman (1960), Deeming (1960), and Griffin (1960). The 
‘‘ Hy-ratios ’’ obtained may be calibrated directly with the aid of Petrie’s original 
standard stars (1953, p. 259). Details of this work will be reported in a later 
paper. 

In conclusion, I would like to say again how much I appreciate all the kind 
help I received,—in the way of observational material, valuable advice, and 
stimulating discussions,—in particular from Professor R.O. Redman and Dr A. D. 
Thackeray, and from Dr R. M. Petrie, Dr R. H. Stoy and Dr K. O. Wright. 


The Observatorics, 
Madingley Road, 
Cambridge: 


1961 May. 
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OBSERVATIONS OF THE SOLAR CORONA 
OUT TO 100 SOLAR RADII 


O. B. Slee 


(Received 1961 June 19) 


Summary 


Observations are described of the scattering suffered by 13 discrete sources 
of radio emission, having ecliptic latitudes between +15° and — 29°, as they 
passed closest to the Sun in the 5-month interval June to October, 1960. 
The equipment used was an E—W interferometer with a base-line of length 
2920A, operating at 85:5 Mc/s. 

It is shown that the average shape of the scattering corona over this time 
interval was elliptical, being more extended in equatorial than polar regions. 
On the average, scattering could be detected out to angular separations of 
about 55 RO and 40 RO in equatorial and polar regions respectively; in 
addition, in the equatorial regions, sporadic increases in scattering occurred 
out to angular separations of 120 RO. Evidence is presented to show that 
the greater equatorial extension of the scattering corona is probably related to 
the active belt on the visible Sun. 

Detailed measurements of coronal scattering obtained by daily observa- 
tion of a relatively intense radio source (MSH 12+08) indicated that signi- 
ficant changes in scattering of a random kind occurred in the interval of 24 
hours between observations. ‘This can be interpreted as due to the relative 
movement of the radio source and a quasi-static system of scattering regions, 
such as coronal rays, if temporal changes in scattering properties can be 
neglected. 

By combining the present results with earlier measurements of the 
scattering angles at closer angular separations, it is shown that the average 
scattering varies as R-*** out to an angular separation of 100 RO in the 
equatorial regions during sunspot maximum. A similar variation in the r.m.s. 
electron densities of the scattering irregularities is suggested if the scattering 
regions can be identified with extensions of the coronal rays. 





1. Introduction.—Up to the present the solar corona has been traced out to 
angular distances of the order of 30-40 R© by means of measurements of the scat- 
tering of radio waves from the intense discrete radio source, ‘Taurus A (for example, 
Hewish (1958), Erickson (1959), Vitkevitch (1960), Hogbom (1960), Gorgolewski 
and Hewish (1960)). However, these observations have been restricted to 
defining the scattering produced in a very limited volume of the corona 
since the apparent path of the radio source with respect to the corona 
changes little from year to year. By using very long base line interfero- 
meters consisting of elements having a large collecting area, which are now 
coming into use, it is possible not only to detect the smaller scattering produced by 
coronal electron irregularities at much larger angular separations but also to 
study its effects on many other weaker radio sources. ‘The latter desirable 
characteristic is due to the fact that the radio Sun, subtending an angle much 
greater than the interferometer fringe spacing, does not produce an interference 
pattern which would obscure or confuse the patterns due to the weaker radio 
sources. The present communication describes observations of 13 discrete 
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radio sources made with an E—W interferometer with a 2920A base line operating 
at 85-5 Mc/s during the months June-November 1960. Since these sources 
possessed ecliptic latitudes ranging from 0°5 to 15° and —o°-5 to — 29°, it was 
possible by systematic examination of their fringe visibilities as they passed 
closest to the Sun, to obtain a fairly detailed picture of the spatial extent of the 
scattering corona. 

2. Observations.—The radio sources used in the experiment are set out in 
Table 1, in which they are identified according to the catalogues of Mills, Slee and 
Hill (1958), Edge et al. (1959), and given I.A.U. numbers where applicable. In 
addition to the necessity of possessing a low value of ecliptic latitude, the sources 
selected needed to exhibit a sufficiently large amplitude 2920A Fourier component 
to make their detection certain above the noise fluctuation level when well clear 
of the solar corona. For the present equipment this requirement was fulfilled for 
sources having a 2920A Fourier component greater than 25 x 10-*° w.m~? (c/s)—1. 
The brightness distributions of the sources are not important in determining 
the scattering introduced by irregularities in the intervening medium; as 
far as the effect of the scattering on the interferometer is concerned, the radio 
source may be regarded as a point having an intensity equal to the amplitude of the 
particular Fourier component. The phase-switched interferometer consisted of 
the E-W arm of the 85-5 Mc/s cross situated at Fleurs, near Sydney, N.S.W., anda 
smaller array of dipoles at a location 1okm due West, which was connected to 


TaBLe I 


Particulars of radio sources observed during the scattering experiments 


Source designation Ecliptic co-ordinates (1950) Observational period 
MSH 3C.—CtCé‘CdC ALL. r B (1960) 
<e 144 o5sN2A 84 ~ sh May 30-July 3. 
06 — 04 161 97 —29 June 7-June 29. 
08 —05 ae 126 —22 June 27—August 2. 
09 +08 230 150 —12°5 July 22—August 29. 
10+02 238 151 — 4°5 July 27—September 12. 
10—02I 249 167 - 7 August 6—September 30. 
11-18 se 181 —15 August 12—October 5. 
12+08 273 185 + 5 August 26—October 31. 
aa 274 12N1A 181 +15 September 4—October 4. 
I2—0l3 275 190°5 — 05 August 31—October 31. 
I2—020 279 : 194°5 + O'5 August 31—October 31. 
13—23 283 204°5 —13°5 September 18—October 31. 
13 —O1I ie 204°5 + 3°5 October 5—November 2. 


the base site via a radio link (Goddard, Watkinson and Mills 1960). ‘The primary 
directivity of the system was 0°-8 arc E-W and 3°-7 N-S to half power points, the 
beam being sufficiently narrow to exclude the intense solar radiation on all except 
the days of closest approach of asource tothe Sun. Pointing of the primary beam 
to the required declination was affected by phasing the dipoles at the remote site ; 
on any one day it was usual to record up to 6 discrete sources over a period of 
about 4 hours around midday, the aerial beam being changed to the appropriate 
declination settings between observations. 

Fig. 1 shows reproductions of some of the records obtained from the relatively 


intense discrete source MSH 12 + 08 while it was in the range of angular separations. 
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19gRO to 109RO. It is clear from these records that for a source giving rise to 
such reasonable signal-to-noise ratios it is possible to measure in some detail the 
variations in apparent size of the radio source from day to day. It should be 
noted that, although there was a general decrease in fringe visibility as the angular 
separations decreased, there were, on many days, large departures from the average 
degree of scattering. Sporadic large increases in the scattering first became notice- 
able when the angular separation was as much as 100 R© and at separations less 
than 60 RO the effects of scattering could be detected on every record. 


225 

















ee | 
| 

































































































































































fit 

















pete 
t 


























+ 


po} = 


Fic. 1.—Reproductions of interference fringes obtained from source 12+08. The angular 
separation for each record is noted near the lower right and generally decreases from left to right 
and top to bottom. The vertical line at the lower left of each record occurs at 124 30™ sidereal time. 
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For the majority of sources used in this experiment, the signal-to-noise ratios 
were between } and } of those shown for the records of Fig. 1; under these con- 
ditions of recording it was possible to say no more than that the source appeared 
about normal intensity or that it was very weak if present at all. It has been 
estimated for such cases that the interference pattern would be classed as ‘‘ not 
present ’’ when the fringe amplitude dropped to below half its normal value. 
Although the weaker sources cannot be used to define the details of scattering in 
the corona, they are still useful as indicators of the average degree of scattering, 
since the points at which their interferometer patterns first disappear and later 
reappear may be used to trace around the Sun a line of constant E-W scattering. 

3. Results—The visibility of the interference fringes for source 12+08 is 
plotted as a function of angular separation from the Sunin Fig.2. Ontheassump- 
tion that the scattered brightness distributions are Gaussian with circular sym- 
metry, a scale has been constructed on the right of the figure to show the corres- 
ponding angular sizes to half-brightness points. We use these equivalent 
angular sizes as a measure of the degree of scattering. 
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The vertical bars inserted at two locations on the diagram indicate the expected 
chance deviations from the curve of best fit and were obtained by computing the 
standard deviations in the amplitudes of fringe patterns due to radio sources far 
removed from the Sun; the half-lengths of the bars represent these standard 
deviations. An examination of Fig. 2 reveals two interesting effects. First, 
departures from the un-occulted mean fringe amplitude became significant on 
sporadic occasions at angular distances as great as 98 RO; the fringe visibility 
was always noticeably reduced at angular separations less than 60RO©. If the 
reduction in fringe visibility is interpreted as due to scattering in terms of the model 
outlined earlier, then it can be stated from these results that the average width of the 
scattered brightness distributions was about 12” of arc to half-brightness points at 
an angular separation of 65 R© and increased rapidly to about 58” of arc at 25 RO. 
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Fic. 2.—Daily measurements of the fringe visibility for source 12+-08 as a function of angular 
separation from the Sun. The full and open circles refer respectively to observations before and 
after conjunction. Tite scale on the right gives the equivalent angular size to half power points. 
The vertical bars inserted at fringe amplitudes of 0-5 and 1-0 have half lengths corresponding to the 
expected r.m.s. chance deviations from the average fringe amplitude. 


No measurements were possible at smaller angular separations due to the loss of 
the interference pattern in the receiver noise fluctuation level. Secondly, we 
note the large dispersion of the points about the average scattering values defined 
by the line of best fit. "These random variations depart from the mean values by, in 
many cases, more than 2 or 3 standard deviations and hence must be regarded as 
indicative of real changes in the properties of the coronal irregularities in the 24 
hours from one observation to the next. ‘The question of whether they should be 
attributed to spatial or temporal coronal variations will be left for later discussion. 

The results for all 13 radio sources observed between the end of May and 
beginning of November 1960 have been plotted in Fig. 3; here, the open circles 
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refer to observations on days when the interference fringes appeared to be 
unaffected by scattering, while the filled circles define days for which the fringe 
amplitude was greatly reduced. The position of the Sun is marked by across and 
the concentric circles are drawn at intervalsof 20 RQ. Over this period the tilt of 
the solar axis with respect to the terrestrial north-south line varied considerably 
and in Fig. 3 the co-ordinate system is fixed with respect to the solar axis. A curve 
has been drawn to enclose the majority of the filled circles and is, as far as can be 
judged, the contour within which the average scattered distribution exceeded 30” 
arc to half power points in the E-W direction. 
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Fic. 3.—Coronal diagram showing scattering suffered by 13 radio sources listed in Table I. 
Open circles refer to observations which show no detectable scattering, while the filled circles denote 
fringe visibilities of about 0-5 or less. The distance scale is defined by circles, concentric with the 
Sun (marked +), having radii increasing in increments of 20 RO. The closed contour, shown dashed, 


marks the average extent of the scattering corona. 

An examination of the coronal diagram shows that the average shape of the 
scattering corona over the 5-month interval was elliptical, being more extended in 
the equatorial regions. In the east equatorial sector, the 30” contour extended 
out to 60 RO, while in the west it appeared to reach about 50 RQ. In the south 
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polar region, the line of 30” scattering almost certainly extended to 40 R© and the 
less numerous data to the north would be consistent with a similar extension. 
The observed axial ratio for the constant scattering ellipse of about 55/40 is 
similar to the general form of the visible corona at this intermediate phase of the 
sunspot cycle. 

Fig. 3 also indicates that significant sporadic increases in scattering were 
recorded on 14 occasions at angular separations in the range 60-120 R© between 
position angles of 0° and 20° with respect to the solar equator. However, on no 
occasion was scattering measurable beyond about 50 R© at larger position angles. 
The significance of the greater extension of the scattering corona in the equatorial 
regions will be discussed in the next section. 

4. Discussion.—The coronal scattering diagram (Fig. 3) showed clearly that 
the contour of constant E-W scattering extended further in equatorial than in 
polar regions; in addition, 14 cases of significant scattering were observed in 
equatorial regions in the range of angular separations 60-120 RO, although no such 
events were recorded in the polar regions. ‘These results are consistent with the 
hypothesis that the electron irregularities which cause the scattering are main- 
tained by the active equatorial belt on the visible solar disk, probably through the 
ejection of material from disturbed regions. 

This suggestion has been tested further by examining the data on the solar 
activity (CRPL, Part B, Solar-Geophysical Data) for the 14 scattering events 
mentioned above and also for 7 cases of abnormally low scattering observed during 
the 5-month interval. Since it was possible that these occurrences were related 
to the presence or absence of active areas near the appropriate limb of the Sun, 
various indications of solar activity were examined over the period concerned. 
There was no significant relationship between the degree of scattering and 
occurrences of important solar flares, intense solar radio emission or severe geo- 
magnetic disturbances ; however, there appeared to be a connection between the 
presence of calcium plages and abnormal scattering occurrences. In order to 
investigate this relationship, a daily limb activity index was constructed consisting 
of the sum of the products of calcium plage area and intensity for plages within 40° 
of each limb. ‘Table II summarises the relationship between the 21 scattering 
events and the computed limb activity index and is indicative of a close correlation 
between them. It should be noticed, however, that this is a one way correlation 
in that the presence of a high activity index was not necessarily associated with 
excessive scattering of the radiation from any particular radio source. Such a 
result is to be expected if each active area on the Sun influences the scattering 


TABLE II 
Correlation between scattering events and solar limb plage activity 
, Above Below 
Scattering 

average average 

Limb activity 

index No. of occasions 
Above average II I 
Below average 3 6 


properties of only a limited volume of the corona so that it becomes a matter of 
chance whether, during any particular observation, the radiation is traversing a 
disturbed region of the outer corona. The latter conclusion is supported by the 
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he detailed observations of fringe visibility for source 12 +08 (Fig. 2) in which the 
mn. scattering angles show a large dispersion about the average for any particular 
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Nebula 1957-1959; A: 26°3 Me/s results of Erickson in 1959; @: 85:5 Mc/s interferometer 
results from source MSH 12-++08 during 1960. 

f It is not possible to decide from these observations whether the changes in 
. scattering in the 24-hour interval between successive observations of the same 
y . . . . . 
i radio source are due to the relative movement of the radio source and a quasi-static 
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system of scattering irregularities or whether dynamic processes related to time 
variations in solar activity play an important part. If temporal variations in the 
scattering properties of the corona over the 24-hour interval can be neglected, it is 
obvious that the physical characteristics of the irregularities in the range 25-120 RO 
vary in an irregular manner over a distance less than 3-5 R© which is the average 
daily change in angular separation. Such a figure would be consistent with the 
transverse sizes of the visible coronal rays if extended outwards to the distances 
considered here. However, in order to account for the observed multiple 
scattering, which depends upon the simultaneous reception of incoherent scattered 
waves from a number of irregularities lying in the plane perpendicular to the line 
of sight, the presence is required of a fine structure having dimensions much less 
than 3:5 RO or 2°5 x 10°km. 

Finally, it is interesting to inquire whether the present scattering measurements 
form a logical extension of the scattering angles already observed at smaller angular 
separations. Measurements of the E—W scattering angles observed between 
5 and 100 RO have been plotted in Fig. 4. Those measurements between 5 and 
15 R© were made by the author by means of a 36’ arc E-W fan beam operating on 
85:5 Mc/s with the Crab Nebula as source of radiation during the conjunctions of 
1957, 1958, 1959; the points plotted between 19 and 100 RO are the present 
measurements on source 12+08. Scattering angles observed by Erickson 
(1959) at 26-3 Mc/s during the 1959 conjunction of the Crab Nebula have been 
scaled to 85-5 Mc/s on the assumption that the scattering angle varies as \? and 
plotted in the diagram to complete the region between 10 and 22R©®. It can 
be seen that the present measurements are entirely consistent with an extrapolation 
of the scattering angles observed at much closer angular separations and suggest 
that, for the 4 years near the maximum of solar activity, the average scattering 
varies as R~-*% out to an angular distance of 100 RO. 

Before inquiring into the corresponding electron densities of the scattering 
irregularities, it would be advisable to examine the factors which influence its 
calculation. It can be shown from the geometrical optics scattering theory 
developed by Chandrasekhar (1952) that the r.m.s. scattering angle, ¢, is related 
to the r.m.s. electron density variation, AN, above (or below) that of the surround- 
ing medium by: 

gb & AN ni? 

where A is the wavelength of the scattered radiation and m the number of irregu- 
larities along the line of sight. ‘The theory assumes that the irregularities possess, 
on the average, identical properties throughout the scattering layer and that the 
scattered brightness distributions display circular symmetry about the direction 
of propagation (i.e. approximately spherical irregularities). It has been shown, 
however, by Hewish (1958) and Vitkevitch, Panovkin and Sukhovei (1959) 
that the irregularities are probably significantly elongated in the radial direction 
from the Sun so that the scattering would be expected to be greatest in the direction 
perpendicular to their long axes. In this case, provided the scattering region is 
small enough so that the directions of the filaments remain approximately constant, 
the above proportionality would remain valid. If the line of sight traverses 
several scattering layers, each having different properties, then the equation must be 
replaced by the appropriate summation. 

It is obvious that in order to estimate r.m.s. electron densities from the observed 
values of the scattering angle it is essential to know the relevant values of n, which, 
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in turn, necessitates a knowledge of the depth of the scattering region, the sizes 
of the irregularities and the spacing betweenthem. Since none of these quantities 
is at present known, it seems futile to attempt calculations of electron density. 
However, if we are justified in identifying the scattering regions as extensions of 
the visible coronal rays, it becomes possible to predict the relationship between 
¢ and AN and how the latter varies, on the average, with angular distance from the 
Sun. 

It seems reasonable to assume that the transverse dimensions of both the ray as 
a whole and its fine structure would continue to increase approximately in pro- 
portion to the radial distance from the Sun. Under these conditions, the value of 
n in the equations would tend to remain constant with the result that AN oc ¢. 
Using the present result that ¢ oc R-** in the equatorial regions during sunspot 
maximum, it seems reasonable to suppose that AN exhibits a similar variation. 


Acknowledgment.—The author wishes to express his gratitude to Mr C. F. 
Fryar whose assistance made these observations possible. 
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RADIAL VELOCITIES OF FUNDAMENTAL SOUTHERN 
STARS, II 


William Buscombe and Pamela M. Morris 


(Received 1961 June 19) 


Summary 


Newly determined radial velocities are given for 15 standard velocity 
stars and 115 other bright stars south of the equator. Measurements from 
two-prism spectrograms (dispersion 90 A/mm at Hy) taken at the Newtonian 
focus of the 74-inch reflector are now entirely consistent with velocities deter- 
mined elsewhere. Most of the present data are from the 30-inch Cassegrain 
spectrograph (dispersion 36 A/mm at Hy). 

Nearly all the stars in this list are in the FK3 and N30 position catalogues, 
of the fifth or sixth visual magnitude. At least three spectrograms of each 
star have been measured. 34 new constant velocities are announced, and 
variability has been detected for the first time in 23 other stars. Observations 
have been resumed on 33 stars previously known to have variable velocity. 
A luminosity class on the revised Yerkes system is assigned where possible. 





This paper continues, in the same form, the data presented for stars observed 
earlier on the programme in Paper I (1). The object of the work has been to 
obtain more reliable radial velocities for certain southern stars whose positions and 
proper motions are known with high accuracy. Of the 198 stars in the FK3 
for which no radial velocity was known in 1940 (2), 133 are south of the equator. 
At least three observations have been completed for each of 34 of these stars at 
Mt Stromlo, and velocities have been published for 21 more by other observatories. 
Progress with the stars having uncertain (usually variable) velocities is slower, but 
orbits have appeared for 7 Mus, ¢ Cen and « Oct. 

A major improvement in the quality of observations has been noticed since 
May 1960, when the faulty temporary camera objective of the three-prism Casse- 
grain spectrograph was replaced. Craftsmen of the Defence Standards Labora- 
tory in Melbourne prepared the new lens system which Mr T. Feuerriegel 
designed specifically for this purpose. 

Although the Zeiss Newtonian spectrograph (3), with its two light flint 
prisms and the 137°5mm camera, gives a dispersion of go A/mm, the optical 
resolution compares favourably with the slower Cassegrain system when the 
optimum spectral purity is maintained behind a slit whose projected image on the 
photographic emulsion is not wider than 18u. Difficulties in keeping this instru- 
ment at constant temperature, experienced during its initial period of operation, 
have been overcome. Rubber insulation on the base casting inhibits heat loss 
by conduction, and a circulating fan and hose prevent thermal stratification of the 
air within the double enclosure. Dr Heinz Gollnow has taken expert care of all 
these instrumental problems. 
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Since the sensitivity of Eastman Kodak emulsions, especially 103aO, available 
during the past four years has declined steadily by more than a factor of two, we 
have been obliged to abandon observation of K stars of sixth magnitude, which 
would now require prohibitively long exposures, and to fill up gaps in the nightly 
observing lists with A and B stars. When a sufficient array of MK standard 
spectra has been observed with the new Cassegrain camera, luminosity classifica- 
tion from the recent plates will be resumed. Plans are in hand to extend this 
work to the G and K stars. 

The practice of observing an internationally accepted star of standard velocity, 
at least once per night when practicable, has been continued. Results of measures 
are presented in Table I not only for Cassegrain spectrograms, which are homo- 
geneous with most of the observations of this programme, but also for Newtonian 
plates photographed with the 74-inch reflector, all reduced with the same wave- 
length system. From the list of lines used in the earlier Cassegrain spectrograms 
(1) only minor omissions were required to obtain consistent results among Miss 
Morris’s measures of the Newtonian spectrograms. In the past two years, no 
deviation from standard velocity, in excess of about twice the internal probable 
error of the mean velocity for a single plate, has been found for either instrument. 
Combined with recent results from the Cape observers, these data may help the 
International Astronomical Union to decide on the suitability of a few additional 
standards considered necessary for the southern hemisphere. 

Following the lists of stars with constant velocities (‘Table I1) and with variable 
velocities (‘l'able III), the individual plate velocities of the double-line binary 
a Oct (HD 199532) are presented in Table 1V. Our preliminary spectrographic 
orbit for this star has recently appeared (4). 


TABLE I 


Standard velocity stars 


HD Star m(HP,,) Spec. V p.e. Obs. Prev. vel. References 
km/s + 
2151 B Hyi 2°80 GrIV +216 o5 17C +21°7a6C Cape Ann. 10 
+23°0 04 9N +23-0a9L Lick Pub. 16 
+23°3a6C M.N. 117, 534 
18884 a Cet 2°54 M2III -27:1 07 7C —258a27 I.A.U.Tr.9 
' —24°4 09 6©2N 
20794 87 Eri 4°15 G5 V +864 ro 5N +88:3b3W Ap. F. 57, 149 
+86:8b6L Lick Pub. 16 
+868 a6C Cape Ann. 10 
+88-3a7St M.N. 118, 609 
+87°5a19C M.N. 119, 638 
29139 a Tau o77 Ks Ill +532 11 4C +54:1a171 IAU. Tr.@ 
62509 6 Gem 1°13 Kol Il + 43 06 3C + 3340167 I.A.U.Tr.g 
80170 120G Vel 5:16 Ks5 — og 04 7N o0a6 LAU. Tr.8 
— 075 a27C M.N. 119, 638 
94510 246G Car 3:78 Ko + 8407 7C + 84aqL Lick Pub. 16 
+7703 7N + 89b3C Cape Ann. 10 
128898 we Cir 3:34 FoV(Sr) + 7:73 06 6C + 71b17L Lick Pub. 16 
+ 76 20 3N + 76a6C Cape Ann. 10 
+ 65a5C M.N. 117, 534 
157457 «x Ara 5:09 KollIl +160 o5 «6©4C) 0+17°6a6 I.A.U. Tr.8 
+16:9a44C M.N. 119, 638 


ee 





XUM 


170’ 


196 


203 


203 


218 


223 


ae ae ee ae ee a 


. 123 
lable 


2, We 
thich 
chtly 
dard 
ifica- 
this 


city, 


sures 
ymo- 
nian 
rave- 
rams 
Miss 
» No 
able 
ent. 
» the 
onal 


able 


nary 
phic 


1G 
8. 








No. 3, 1961 Radial velocities of fundamental southern stars 
TABLE I—continued 
Standard velocity stars 
HD Star m(HP,) Spec. VV p.e. Obs. 
km/s + 
171759 ¢ Pav 3°94 Ko —16°3 o8 5C 
—15°7 04 O6N 
196171 « Ind 3°05 Kolll — 12 o2 2N 
203608 y Pav 4°22 F6V —29°3 04 4N 
203638 33 Cap 5°27 Kolll +2372 o2 3C 
+2271 o3 2N 
218594 88 Aqr 3°75 Kolll +228 o5 4C 
223647 y! Oct 4°99 G5 Ill +1375 o2 3N 
TABLE II 
Stars with constant velocities 
HD Star m(HP,,) Spec. V p.e. Obs. 
km/s + 
2363 69G Scl 5°77 G5 —~% ops 
8810 9G Hyi 5°62 Ks5 +23 1I°§ 3 
*9362 5 Phe 3°84 Ko — 8 16 9 
10476 107 Psc 518 KiV —34 1°4 2 
11937 x Eri 361 GsIV —- 6 Io 3 
12363 o Hyi 6-15 F2 +6 04 3 
12477 26G Hyi 5°97 Ka + 5 04 4 
12573 201G Cet 5°59 As +13 O72 3 
14728 239G Cet 5°77. Ko — 3 07 3 
15248 x Hyi 5°84 Ko +22 08 3 
15471 29G For 5°95 Gs — 2 O85 3 
16538 v For 5°67 Fs +28 o7 3 
16975 40G For 5°97 Gs +17 16 3 
20640 40G Hor 5°68 Ko — 9 04 3 
23719 51G Hor 5°44 Ko — 4 08 2 
24188 61G Hyi 6°46 Ao + 2 06 3 
25371 163G Eri 5°78 AoV +21 13 4 
27256 a Ret 3°26 GélIl +35 04 3 
27563 210G Eri 577 Bslll +8 13 5 
29116 v Men 5°72 Fo +14 19 4 
30127 264G Eri 546 ArV +5 144 
30432 15G Cae 5°82 Ko -— 6 13 3 
32515 26G Cae 5°78 Ko +29 1'O 3 
33256 68 Eri 5:10 Fs V +7 FO 2 





wo 


Prev. vel. 


—16-9b8L 
—16-9b3C 
—18agL 
oo b6C 
— o8a5St 
—301a10L 
—31°5b3C 
—28:8a5C 
—30°3.a3St 
+22:0a7 
+22:1a39C 
+2r1a6L 
+21°4a5C 
+14°5b3L 
+12°0a4St 
+13:93a9C 


Prev. vel. 


—65a8L 
-—-71b7C 
—34b4W 
—33a3L 
—6c8L 
—6a15C 


+g9d1D 


-—2di1L 


+36b7L 
+36b9C 


+10a3L 
+ 8b3C 
+10a2W 
+24°7a19 
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HD 


40972 
*45348 
48915 


*50877 
54605 
56855 

*57623 
59612 


68434 
69267 


*74006 
89388 
90853 


119938 
123123 


139233 
141891 


154783 
159877 
163917 


165259 
182509 
188584 
197725 
199684 
200702 
201647 
203532 
209529 
210853 
214632 
221323 
224930 
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80G 
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187G 
196G 
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10G 
5G 
85 
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Car 
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Vol 
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Car 
Cnc 


Pyx 
Car 
Car 


Cen 
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Ser 
Oph 
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Tel 
Pav 
Cap 
Mic 
Mic 
Mic 
Oct 
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Oct 
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TaBLE I[I—continued 


m(HP,,) Spec. V p.e. Obs. Prev. vel. 
km/s + 


5:90 Ar V +12 06 4 
—o'79 Folab +18 04 7 =+20°5a114 


—1'68 ArV — 6f 06 3 — 8ar1orL 
3°79 K3lab +37 03 3 +37b5L 
+36b4C 
1'90 F8 Ia +34 12 4 +35a46L 
+33a20C 
244 KollII +16 06 2 +15:5a10L 
+16-1a7C 
3°94 F8II +22 18 6 +22°6a6L 
+21:7b3C 
4°74 A3 +40 2:00 § +37c8L 
+36c3Y 
5°75 Az +25 O07 3 
3°42 Ka +23 1°5§ § +22b5L 
+23b3C 
+19a2V 
+18 4W 
3°87 G5 —12 0613 —148a9L 
—15°4b3C 
3:23 Ksilb + 6 062 + 8agL 
+10b3C 
400 Foll +10 06 5 +10birL 
+ 6a3C 
6:01 A3 +16 o8 4 
3°22 K2III +27 o7 2 +27a19L 
+28b6C 
6°56 Bo +5§ 25 
2°97 F2V Oo 83 7 oai5L 
+ 3b6C 
5°83 A3 —30 O05 4 
5°83 As —- 8 10 4 
3:22 GolIII +14 04 2 +13a7L 
+13a6C 
5°84 Fs +1I O'2 3 
5°39 Kz ~s oy 9 
5°57 Ko +28 o7 4 
' 5°86 Ao +18 04 3 
6:05 Fo +15 04 3 
5°34 Ko +1II 04 3 
5°76 Fs +12 12 2 +11d1W 
6°37 Bz IV © O77 4 
5°40 Ks —II 12 4 
5°58 A6GIV +17 Ir 3 
5°75 Ko —32 06 3 
5°70 Ko + 8 o5 4 
5°72 G2V —32 09 2 —33¢5W 
—35a3V 
—37 SY 


* Possibly slightly variable in velocity 
t+ Epoch 1954.1. 
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Prev. vel. 


oc6L 


— 3d4L 


242 
TaBLe [11—continued 
HD m(HP ») JD V p.e. 
Star Sp. 7. km/s + 
224990 4°91 4998: 106 +3 49 
£Scl Bs III 7173°170 + 9 3°6 
7211°048 +21 3°9 
7218128 +13 5°2 
7245°967 +8 48 
225253 5°54 5361°123 +12 FF 
45G Tuc B8 III 5380°081 - 5 2°8 
6036°251 + 3 3°5 
* Orbit. 
TaBLe IV 
Radial velocities of HD 199532 
JD V; Ve JD 
7 re 
5314°129 +15 6353°222 
5356°905 +25 6353°336 
5745°997 +10 6379°181 
6036°134 + 7 6393°178 
6038-118 +86 6406°151 
6067-068 +66 6410°122 
6084:013 +92 — 2 6711°307 
6084:°088 +88 — 4 6711°317 
6085-969 + 34 6718-283 
6086-087 +60 6720°267 
6092°025 +64 +20 6737°231 
6120°923 +82 6778-137 
6154°920 +32 68037045 
6182937 +66 
TABLE V 
Observers and measurers 
Observed 
L. Aller* 5 
W. Buscombe 366 
H. Gollnow 67 
J. Graham ° 
W. Heintzt 8 
P. Morris 27 
A. Przybylski 84 
D. Sher 19 
J. Whiteoak 5 
Totals 581 





* Michigan University: 
+ Munich Observatory: guest observer, 1954-55. 


visiting professor, 
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A STUDY OF THE GALACTIC CLUSTER IC 2602 
I. A PHOTOELECTRIC AND SPECTROSCOPIC INVESTIGATION 
- John B. Whiteoak 


(Communicated by A. W. Rodgers) 


(Received 1961 June 20) 


Summary 


The results are presented of a photoelectric and spectrographic study of the 
galactic cluster IC 2602. In the investigation, observations of 29 probable and 
40 possible members yield a distance of 155 parsecs and cluster age of 10” years. 
It is a poorly-populated cluster and its luminosity function shows a deficiency 
in the number of stars of spectral type F compared with that predicted by the 
Initial Luminosity Function. It is suggested that the star-group is a local 
condensation of the Scorpio-Centaurus association. 





Introduction 


The galactic cluster IC 2602 (Melotte 102) is an open group of stars distributed 
about the Bo star Theta Carina. The brightest members are in a group 13° in 
diameter centred at «=10"42™, 5= —64°(1950) (J4'=289°7, bUL= —4-°9). 
The cluster is shown in Plate 7: a 3-min exposure in visual light, taken with the 
20/26-inch Schmidt telescope of the Uppsala Southern Station. The stars with 
assigned numbers are those measured in the search for cluster members. The 
small dense group of stars on the outskirts of IC 2602 at «= 10" 41™, 5= — 64°-9 
is the more distant cluster Melotte 101. 

In a recent photoelectric study of the cluster, Braes (1) derived a distance of 
150 parsecs, a value less than those of former investigators (2). From the absolute 
magnitude of the brightest main sequence star Theta Carina, and the position on 
the main sequence that gravitational contraction occurs, he derived an age of 
15 x 10° years. 


Observations 


1. Photoelectric.—Three telescopes were used to obtain the photoelectric 
observations of the suspected cluster members. The investigation was begun 
with the Mt Stromlo 20-inch f/18 reflector, but most of the observing was carried 
out on the Mt Bingar 26-inch f/12 reflector and the Mt Stromlo 50-inch f/18 
reflector. The photoelectric equipment consisted of an EMI photomultiplier, 
replaced by a refrigerated RCA 1P21 for the more recent observations, a GR 
amplifier or one built at Mt Stromlo, and a Brown Recorder. The filters used 
were those recommended by Johnson (3) for the linear transformations from 
instrumental magnitudes to those of the UBV system. The combinations of 
filters and photomultipliers are shown in the following table. 
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Telescope Cell Filters 
20-inch 1P21 OY4; BG12+GG13; UG2 
26-inch EMI; 1P21 C3384; C5030+ GG13; C9863 
50-inch 1P21 C3384; BG12+ GG13; UG2 


C3384; C5030+ GG13; C9863 


As mentioned above, the initial observations were carried out on the 20-inch 
telescope, but rarely were the nights of sufficient quality to compare accurately 
the cluster stars with the equatorial UBV standards of Johnson and Harris (4). 
Magnitudes and colours measured with this telescope were given lower weight 
than those obtained elsewhere. 

The 26-inch telescope was used primarily for the transformation of the UBV 
system from the equatorial standards of Johnson and Harris to several stars in the 
vicinity of the cluster. A total of 23 UBV standards was measured for the 
transformations, the average number observed per night being 12. From five 
nights of good quality, the transformation equations determined were : 


V —V.,= —0°102(B—V)+0°142 


B-V 1°156(B— V )og + 1°964 
U-B 1-051(U— B)s,— 1-051 


with a mean error per measurement of o-o12 in V, o-o15 in B—V, and o-o18 
in U—B. 

The sensitivity of the photomultiplier was checked during each night with a 
radium source fitted in the filter wheel. To obtain extinction coefficients each 
night as well as to check on the constancy of transparency throughout it, two stars, 
one of early, the other of late spectral type, were observed several times. For 
nights of doubtful quality, the following mean extinction coefficients were used : 


K, =0°'18 
Kp_y =0°15-0°04 (B—V) 
Ky-,=028 


In general, stars were observed through focal apertures of 16” and 29”. 

The observational procedure with the 50-inch was the same as with the 26-inch, 
except that stars were observed with respect to the secondary standards. The 
following mean extinction coefficients were used on night of doubtful quality : 


Ky =0O° 18 
Kp_y =0°12-0°04 (B-—V) 
Ky -p=031 


The reduction of observations was effected in the manner recommended by 
Johnson and Morgan (3). After the installation of an IBM 610 computer at 
Mt Stromlo early in 1960, subsequent reductions were carried out on it. 


2. Spectrographic.—To obtain a luminosity and spectral classification of some 
of the cluster stars, spectrograms taken by Drs H. Gollnow and A. Przybylski 
for the measurement of radial velocities were kindly made available. ‘They were 
obtained with a Zeiss two-prism spectrograph mounted at the f/5 Newtonian 
focus of the 74-inch reflector. The dispersion of the spectra is go A/mm at 
H-gamma. With IIaO and 103aO emulsions, and a projected slit-width of 20 
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microns, exposure times varied from a few seconds for the star Theta Carina to 
100 minutes for a ninth magnitude F2 star. 

The spectral and luminosity classification of the stars was effected by compari- 
son of their spectra with those of some of the MK standards listed by Buscombe 
(5) that had been obtained with the same spectrograph. It was based on the 
ratios of spectral-line intensities shown in the MKK Atlas (6) to be sensitive to 
spectral type or luminosity. 


3. Photographic.—In a search for fainter possible members of the cluster, 
plates centred on a=10"41™, 6= —64°:2 (1950) were taken in blue and visual 
light with the f/3-5, 20/26-inch Uppsala Schmidt by Dr B. Westerlund and Mr D. 
Sher and kindly made available to the author. Plates in each colour were given 
exposures of 20 seconds, 30 seconds, and 3 minutes. The limiting stellar magni- 
tude on plates with the last exposure exceeds 17. The photographic emulsion- 
filter combinations were : 


Visual plates 103aD + 2mm GG14 filter 
Blue plates IlaO +2 mm GG73 filter 


Images of selected stars on the plates were measured with a Haffner iris 
diaphragm photometer. Only one setting per star per plate was necessary as the 
setting errors are less than those due to the plate emulsion. The photometer 
readings were converted to UBV magnitudes on the basis of those stars already 
measured photoelectrically. To correct for the large field errors occurring on the 
plates, a basic calibration curve was constructed with the stars of Mel. 101 observed 
with the 50-inch, and the field errors determined from previously-measured stars 
that are on different parts of the plate. The mean error of a single reduction, 
o™-05, although large, is sufficient to determine whether the star is a possible 
cluster member. 

Approximate spectral types of the stars in and near the cluster were shown on a 
plate of the field taken with the Schmidt telescope fitted with a 4° objective prism. 
With an exposure time of 10 minutes, spectra with a dispersion of 480 A/mm at 
H-gamma of stars as faint as 12th magnitude could be classified. An approximate 
spectral classification could be achieved by the comparison with the spectra of stars 
classified in the Henry Draper Extension (7). 


Results 


The stars measured photoelectrically are listed in Table I. The columns 
are: (1) the star’s number in the present study; (2) its HD or HDE designation ; 
(3) right ascension for the epoch 1950; (4) declination for the same epoch; (5) V; 
(6) B—V; (7) U-B; and (8) the number of times that the star was observed. 

The stars having a great number of observations were those used as secondary 
standards. 

The stars 1-33 were compiled previously at Mt Stromlo in a list of probable 
members for a radial velocity study of the cluster. The remaining stars are some 
of comparable brightness situated in the same general direction as the first group. 
Band V measurements of 27 of the listed stars have been made by Braes. The mean 
differences are (Braes-Whiteoak) : 


in V —o™917 (+0°018) 
in B-—V +0™-013 (+0-014). 








248 John B. Whiteoak Vol. 123 No. 
TaBLeE I 
Stars observed photoelectrically 
N 
1950 4 
No. HD/HDE R.A. Dec. V B-V U-B n 5 
I 92385 104 36-2 — 64°47 6°75 —0:08 — 0°30 7 5: 
2 92467 10 37°1 —64 14 6°99 0°03 —o'16 8 5 
3 92478 10 37°3 —64 42 7°56 0°05 0°00 9 5: 
4 92492 10 37°4 —63 25 9°27 or12 0°04 II 5! 
5 92535 10 37°6 —63 32 8:25 o'21 o'10 9 5 
6 92536 10 37°6 —63 51 6°34 —o0'08 —o°31 8 6 
5 92568 10 37°9 —63 26 8-58 0°42 o'l4 9 6 
8 92569 10 37°9 —64 26 9°46 0°22 Orls 6 6 
9 92570 10 37°9 —64 50 8-86 0°47 0°02 I 6 
10 92664 10 38°5 —64 50 5°52 —O'I7 —0'59 7 6 
II 92715 10 38°9 —64 24 6°82 — 0°02 —O'l4 9 6 
12 92783 10 39°3 —64 13 6°74 — 0°05 — 0°22 13 6 
13 92837 10 39°8 —64 00 7°18 0°00 — 0°08 13 
14 92896 IO 40°! —63 I5 732 0°22 0°06 7 
obs 
15 92938 IO 40°5 —64 12 4°81 —O'15s —o'58 20 
16 92966 10 40°6 —64 08 7°28 —o'ol —o'10 14 To 
17 92989 10 40°8 —64 26 7°58 0°04 —orol 8 pla 
18 93030 10 41°2 —64 08 2°78 —0'23 — 1°03 8 U- 
19 93098 IO 41°7 —63 48 7°61 0°04 0'00 9 it a 
20 93163 IO 42°1 —63 59 5°79 — 0°02 —0'55 II 
21 93194 10 42°3 —63 42 4°84 —o'l4 — 062 9 
22 93209 IO 42°4 —64 25 9°41 0°32 0°25 10 
23 93405 IO 43°7 —63 58 9°17 0°44 o-o! II 
24 93424 10 43°8 —64 27 8-14 Orls 0°08 14 
25 93517 IO 44°4 —64 20 7°86 0°09 0°07 12 
26 93540 IO 44°5 —64 15 5°36 —o'lo —o0°48 fe) 
27 93549 IO 44°7 —64 00 5°26 —o'09 —o'48 12 
28 93600 IO 44'9 —64 22 8°44 0°50 0°02 15 
29 93607 10 45°0 —64 07 4°87 — Ors — 0°65 9 
30 93648 10 45°3 —64 00 7°85 o'll 0°08 10 
31 93714 10 45°8 —64 17 6°57 0°00 —o'60 7 
32 93738 10 46:1 —64 00 6°48 0°00 —o'16 9 
33a ©. 93796 10 46°3 —63 21 10°28 0°63 o'13 3 
b 9°94 0-58 0°07 3 
34 91906 10 33°1 —63 53 7°46 0°03 oO'ol 2 
35 307860 10 33°2 —63 52 8:24 0°22 0°03 I 
36 =. 307861 10 33°2 —63 57 10'OI 0°50 0°22 2 
37. 91944 16 33°4 —63 58 8:68 0°05 — 0°38 2 
38 = 91997 10 33°8 —64 02 8-95 0°04 — 0°33 4 
39 92066 IO 34°2 —63 59 8-40 0°00 —0°43 4 
40 307859 IO 35° —63 52 10°35 0°20 0°05 I 
41 307842 10 35°5 —63 24 9°33 0°73 0°33 I 
42 307844 10 36°8 —63 19 9°30 1°41 1°47 I 
43 307856 10 37°7 —63 50 10°58 1°22 I°ls 2 
44 92663 10 38°6 —63 45 7°81 T°S 1°68 17 
45 307926 10 38°6 —63 46 9°70 1°34 1°49 6 
46 307919 10 39°3 —63 30 10°06 I-14 0'90 I 
47 92821 IO 39°7 —63 31 8°87 1°28 1°31 I 
48 307921 IO 40°! —63 38 9°04 1°58 1°60 2 
49 307922 IO 40°4 —63 42 9°99 0°60 o'18 2 
50 307924 IO 40°7 —63 44 10°70 0°62 0°02 I 
51 307931 IO 40°9 —63 56 IO°II £22 1°07 14 : 
52 307937 IO 41°4 —64 o1 9°76 Ort —0'26 4 
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TABLE I—continued 
1950 
No. HD/HDE R.A. Dec. V B-V U-B n 
53 307933 10 41°5 —63 56 8-90 0°43 0°25 I 
n 54 307934 IO 41°7 —63 58 10°08 1°63 1°90 10 
7 55 93115 10 41°8 —63 24 7°85 1°52 1°78 6 
8 56 93269 10 42°8 — 63 33 8-15 I°l4 1°04 6 
9 57 307963 10 43°9 —64 31 10°06 1°15 0°89 3 
rT 58 308006 10 44°6 —63 35 9°78 1°53 1°56 5 
9 59 308015 IO 44°7 —64 05 8°95 I°IQ 0°85 8 
8 60 308012 IO 44°7 — 63 53 10°06 0°49 0°03 2 
9 61 93777 10 4674 —63 38 8-51 1°16 1°04 3 
6 62 93874 10 47°0 —63 35 8-20 O17 0°09 2 
I 63 93892 10 47°1 —63 48 8-99 0°50 0°06 I 
7 64 94066 10 48°1 —64 17 7°90 0°08 —0°38 5 
9 65 94174 10 48°8 — 64 13 7°76 orl 0°06 3 
3 66 94422 10 50°6 — 64 09 8:24 o'ro —0°39 2 
3 Fig. 1 shows the relation between the colours B— V and U—B of the stars 
ef observed. ‘The stars 1-33 are represented by open circles, the remainder by dots. 
4 To the former has been fitted the unreddened colour relation given by Arp (8), dis- 
8 placed in B—V by o™-04, the mean colour excess of the cluster, (Ez_y), and in 
S U-—B by o™-03. From the positions of several of the circles relative to the curve, 
: it appears that the stars whose numbers are 7, 20, 22 and 31 may not be members, 
9 
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Fic. 1.—Observed relation between colours of stars. 
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assuming they represent normal main sequence stars. Of these, Braes has shown 
that 20 and 31 have radial velocities differing considerably from the mean value for 
the cluster. Most of the dots represent either more-reddened main sequence 
stars or background giants of late spectral type. Stars that may be associated 
with the cluster are Nos. 34, 35, 36, 41, 43, 49, 50, 60, 62, 63, 65. From the 
relation Ay =3:0E;_,, the mean visual absorption for the cluster (A,;) is 
o*-12. 

Table II contains the results of the spectral and luminosity classification of 
the stars. For the stars observed spectrographically with types earlier then Ao, 
values of (B— V4.) were derived by the ‘‘Q’’ method (Johnson (g)). For stars 
of later spectral types, Arp’s calibration (8) was used. 

The columns of the table are: (1) Star number; (2) spectral type and lumino- 
sity class; (3) The number of spectra examined; (4) Spectral type from the 
““Q” method; (5) Ep_y; (6) (B—V)y; (7) Vo; and (8) Notes on the spectra 
examined. 

Tas_e II 


Determination of cluster reddening 


No. Sp. n_ Sp. (“Q’’) Ep-y (B-V)po V, Notes on spectra 
I Bo V 4 B8 “00 — 0°08 6°75 
a Bo IV 5 Bo “09 —0-'06 6°72 
3 Ao IV 4 Ao 06 —o'ol 7°38 
4 Ao IV 4 Ao “12 0°00 8-91 non-member. 
5 A3V 4 “12 0°09 7°89 
6 B8 IV 5 B8 ‘Ol —o'09 6°31 
7 A7p 5 0°38 8-46 enhanced lines of Sr 
and Fe. 
8 Bos IV 7 Ao “23 —orol 8-77. non-member. 
9 F7V 4 0°43 8-74. Close visual binary. 
F5IV 5 
10 Bo p 5 Bs —o'2I 5°40 ©Sistar. 
II Bos V 4 Bo'5 "02 —0'04 6°76 
12 B8 V 5 Bo ‘OI —o:06 6°71 
13 Bo V 4 Bo's 03. — 0°03 7°09 
14 As IV 5 ‘O07 ols Fey 
15 B3 V 5 Bs ‘Ol —o'16 4°78 
16 Bo's5 V 2 Bo's5 "02 —0'03 7°22 
17 A1rlV 6 Ao O05 —o-ol 7°43 
18 Bo V 4 Bo 05 —o028 2°63 Og'5 (V). 
19 Ao V .. Ao 05 —o'ol 7°46 
20 B3 V 4 B4 "15 —O'17 5°34 B3:V(V);non- 
member. 
21 Bs V 4 B4 03 —0'17 4°75 Bs Vn(V). 
22 A4 IV 5 20 or12 8-81 non-member. 
23 F2V 4 06 0°38 8-99 
24 A4 IV 5 03 orl2 8:05 
25 A2V 4 03 0°06 Y hag i) 
26 B6 V 5 B6 04 —O'14 5:24 B7:V(V). 
27 B7 IV 4 B6 "O05 —o'14 Ce 8 
28 F7V 4 ‘00 "50 8°44 
29 B4V 6 B4 03 —o'18 4°78 B4IV(V); B3 IV (M). 
30 AzV 4 "05 0°06 7°70 
31 B2 IV 8 B3 "24 —0'24 5°85 non-member. 
32 Ao V 5 Bo ‘05 —0°'05 6°33 
33a G2IVv 2 ‘00 0°63 10°28 
b F8&V 2 05 0°53 9°79 
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In Col. (8), previous classifications are shown. ‘Those followed by (V) 
were observed by Mme de Vaucouleurs (10), those by (M) by Miss Morris (11). 
From Col. (5), one can see that the stars Nos. 8, 20, 22 and 31 are more 
reddened than the remainder and are presumably background stars, as they 
appear to have normal spectra. Stars with peculiar spectra, Nos. 7 and 10, have 
been corrected for absorption by the reduction of their values of V and B—V by 
the means of individually-determined colour excesses and visual absorptions, 
o™-o4 and o™-12 respectively. In addition to the four excessively-reddened stars, 
No. 4 is too faint for its spectral type for it to be a member. The stars observed 
spectrographically are represented in the colour-magnitude array, corrected for 
absorption, of Fig. 2. ‘The doubtful members are shown as open circles. ‘To the 
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Fic. 2.—Colour-magnitude array of the cluster, corrected for absorption. 


filled circles has been fitted the ‘zero-age ’ main sequence of Johnson and Iriarte 
(12), displaced 5™-94 vertically, the distance modulus of the cluster. It is equiva- 
lent to a distance of 155 parsecs, similar to the value (150) determined by Braes. 
A feature of the main sequence is that it does not seem to extend to magnitudes 
fainter than V,=8-7. This could be due to the effects of gravitational contraction 
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for the fainter stars, a lack of stars of spectral type F in the cluster, or the fact 
that the survey for cluster members is incomplete at the fainter magnitudes. 

A spectroscopic parallax for the cluster was not derived, for the MK—M, 
calibration would not account for the fact that the brighter stars have already 
evolved to a position above the main sequence and are therefore more luminous 
than their main sequence calibration would predict. ‘The same argument applies 
to the apparently contracting stars, although until some other membership 
criterion like radial velocities can be considered, it is not certain whether they are 
contracting cluster stars or foreground objects. 

On the visual, blue and objective-prism Schmidt plates, an area of sky 12° in 
diameter was examined in a search for fainter possible members. V and B—V 
measures were obtained for such stars. ‘The additional possible members are 
listed in Table III. The columns are: (1) the number of the star in the present 
survey; (2) the HD or HDE number of the star; (3) V; (4) B—V; and (5) 
its spectral type—either that assigned to it in the Henry Draper Extension or that 
derived from the spectrum on the objective prism plate. 


Taste III 
Possible members observed photographically 


No. HD/HDE V B-V Sp. 
67 10°75 I'lg K 
68 II‘ol 1°44 K-M 
69 308023 9°80 0°96 Ko 
7° 307849 10°43 0°33 F 

71 10°54 0°40 Gs 
72 307845 10°90 0°50 Go 
73 307885 10°83 riz Ko 
714 308011 9°82 1°16 Ko 
75 10°42 0°53 Go 
76 307942 10°03 0°42 Fs 
77 397944 9°40 1°07 Ko 
78 10°93 0°69 Go 
79 11°46 0°83 G 
80 307918 10°72 0°88 G 

81 307938 10°48 0°57 Go 
82 307929 10°65 1°03 K2 
83 10°96 I‘Io G7 
84 307960 11°33 o'9I G5 
85 307936 10°87 0°70 F7 
86 307955 10°84 0°95 Ko 
87 307969 10°98 1°32 K7 
88 307979 10°96 0°68 Go 
89 307956 10°96 0°64 G 
go 308016 10°59 0-98 K2 
gI 10°64 0°56 Go 
92 11°08 o'51 Go 
93 10°87 1°55 Ks 
94 10°12 1°26 K2 
95 10°70 1:08 Gs-Ko 


The total population of the cluster is shown in the colour-magnitude array of 
Fig. 3. The probable members are shown as filled circles, additional photo- 
electrically-measured stars as open circles, and those measured photographically 








as d 


add 


mag 
lum 
sim 
sev 
lack 





» fact 


eady 
nous 
plies 
ship 
y are 


¢° in 
—V 
$ are 
sent 


| (5) 


that 


r of 
ito- 








No. 3, 1961 A study of the galactic cluster IC 2602 253. 


as dots. Too the first group has been fitted the ‘zero-age’ main sequence. The 
additional curve leaving the main sequence at V=8-4 would define the bright 
magnitude limits at the faint end of a colour-magnitude array of a cluster with a 
luminosity function similar to that of NGC 2264 (Walker (13)) but with an age 
similar to that of IC 2602. Despite the fact that in addition to fainter members, 
several field stars probably are included in the figure, the cluster still shows the 
lack of F stars mentioned previously. 
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Fic. 3.—Colour-magnitude array of the cluster in which fainter stars, possibly members, have 
been included. 


Discussion 


Luminosity function—The luminosity function of the cluster is shown in 
Fig. 4, the distribution of the number of stars per unit magnitude (N[My;]) 
with absolute magnitude (M;,)._ The broken line for My greater than 3 represents 
that part of the luminosity function which is not so well-defined. The curve in 
the figure is the Initial Luminosity Function adopted by Limber (14). It was 
normalized to the cluster so that the number of members it predicts between the 
limits of absolute magnitude —1 and +2, a range within which the observed 
membership of the cluster is complete but has not been affected yet by evolution, 
is equal to the observed number. Compared with the Initial Luminosity Func- 
tion, that of the cluster appears to be lacking in stars of absolute magnitude +3 
and +4. This is consistent with the conclusions of van den Bergh and Sher (15) 
that galactic clusters do not have a unique luminosity function. Attempts were 
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made previously to estimate the number of faint cluster members by star counts 
and the distribution of stars of late spectral types, i.e. stars with objective prism 
spectra showing pronounced H and K lines. Both were unsuccessful due to the 
irregular background absorption. 
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Fic. 4.—Luminosity function of IC 2602. 


Age of the cluster.—The age of the cluster can be derived by two methods: 
(1) from the absolute magnitude of the brightest star on or near the main sequence; 
and (2) the position on the main sequence below which the stars are still in the 
process of gravitationally contracting on to it. 

The brightest star near the main sequence, Theta Carina, at present has an 
absolute magnitude of — 3-3. It is already evolving from the main sequence and 
if it is assumed that it has moved almost vertically from it, its absolute magnitude 
on the main sequence would have been — 2-8. According to the results of recent 
calculations of star models (Limber (15)), such a star would attain an absolute 
magnitude of —3-9 in a time of 1-2 x 107 years since its formation before the 
hydrogen in its core js exhausted and it evolves from the neighbourhood of the 
main sequence. As the time that it takes to brighten to this magnitude from its 
value on the main sequence is small compared with the age of the star, then 
assuming that all the cluster stars are coeval, the cluster would have an age of 
1'2x 10’ years. ‘This value is greater than that determined with the older 
Schénberg—Chandrasekhar theory. 

The cluster main sequence appears to terminate at V,=8-7 i.e. Mp= +2°8, 
equivalent to a star of mass 1} solar masses. Henyey, Le Levier and Levee (16) 
have calculated that a star of such a mass contracts gravitationally on to the main 
sequence in a time of 8 x 10% years. If it is assumed that the termination is due 
only to the effects of gravitational contraction, which as mentioned previously is 
doubtful, then the age of the cluster would be 8 x 10° years. For IC 2602, the 
age calculated from this method would be expected to be less accurate than the 
first. 
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From both considerations, it appears that the age of the cluster is about 
10° years. 

Proper motions.—Proper motions of eighteen of the cluster stars are listed by 
Boss (17). ‘The mean values are: 

— 08-0024 + 0:0008 ; —0":005 + 0-010. 
at a distance of 150parsecs in the direction of the cluster, the apparent proper 
motions due to solar motion are: 

— 08-0015 ; +0”-007. 

It appears that most of the apparent motion of the cluster is due to the effects 
ofsolar motion. For this reason it is difficult to detect members using the criterion 
of common proper motions—the General Catalogue lists many stars, obviously 
non-members of the cluster but in its general direction, that have comparable 
proper motions. 

From an analysis of the proper motions of the members, Markarian (18) 
suggested that the cluster is expanding. A similar conclusion in the case of the 
association I Lac was criticized by Woolley and Eggen (19), and it is quite likely 
that for IC 2602, any apparent expansion effect is due also to systematic errors in 
the proper motions, particularly as the determined expansion age (3 x 10° years) 
is considerably smaller than the age derived in the present study. 


The relation of IC 2602 to the Scorpio-Centaurus Association 


As the Scorpio-Centaurus association is a close group of young stars, it is of 
interest to investigate the possibility that IC 2602 is related to it. According to 
Blaauw’s list of certain members (20), the association extends from og" oo™ to 
19"00™ in right ascension and in the direction of the cluster it has a similar 
declination. The mean of the radial velocities of the cluster stars listed 
by Wilson (21) is +24km/sec, a value comparing favourably with that of the 
association in the direction of the cluster (Buscombe and Morris (22)). 

From a study of the proper motions of the association stars, Bertiau (23) 
derived a distance of 170 parsecs. Moss Morris (11), on the other hand, deter- 
mined a distance of 220 parsecs from a spectral and luminosity classification of the 
certain members. She used the relation between MK classification and absolute 
magnitude of Johnson and Iriarte (12), where the luminosity class V is derived 
from the stars in the vinicity of the Sun. For stars of early spectral type, it 
includes those stars that are already beginning to evolve from the main sequence, 
so that the mean absolute magnitude for a star of a certain spectral type is more 
negative than the value it had on the main sequence. ‘The use of sucha calibration 
would result in an overestimation of the distance. A colour-magnitude array 
of the association, constructed from the magnitudes and colours of Bertiau for 
Blaauw’s probable members, reveals that the average reddening (excluding the 
more heavily reddened star Sigma Sco) is o™-04 and the average distance of the 
group is 150 parsecs. ‘The latter value is only approximate because of the lack 
of stars of later spectral types. Both of these figures agree with those of the 
cluster. 

The age of the association, determined from the main sequence stars of earliest 
spectral type (Bo), is about 10’ years, similar to that of the cluster. 

The above evidence suggests that IC 2602 may be a local condensation of the 
Scorpio—Centaurus Association. 
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CAN GRAVITATIONAL FORCES ALONE ACCOUNT FOR GALAXY 
FORMATION IN A STEADY STATE UNIVERSE? 


II. BonpI-GOLD UNIVERSES 


Martin Harwit 
(Received 1961 July 10) 


Summary 


The accretion problem in a Bondi—Gold steady state universe is examined, 
We show that if galaxies are to be formed in a homogeneous universe by 
predominantly gravitational means, the cosmic density must be at least 

Pmin = 15/47yT? 

where y is the gravitational constant and T~! is Hubble’s constant. This is 
ten times the value initially specified by Hoyle for a steady state universe whose 
expansion rate is governed by the general relativistic field equations. For a 
Hubble’s constant 2:5 x 107*8 sec™!, the minimum density pmin=9 x 107° g 
cm-*, This is higher than the usually estimated maximum density of the 
Universe, but present estimates are not yet well enough established to 
rule out definitely gravitational galaxy formation in an arbitrary steady state 
universe. 





1. Introduction.—In 1948 two different types of steady state cosmologies were 
proposed in quick succession. Bondi and Gold (1) described a steady state 
(BGSS) universe which was not necessarily governed by the laws of general 
relativity. ‘The authors felt that there was no reason to believe that general 
relativity was valid on the scale of the entire universe; observations had only 
verified the theory on the scale of the solar system. 

An alternate version of the theory proposed by Hoyle (HSS) (2, 3) incorporates 
the general relativistic field equations and is, therefore, able to relate the expected 
mean density py of the universe to 1/7, Hubble’s constant. The relationship is 

Po = 3/8ryT”. 

A recent (4) examination of gravitational galaxy formation in steady state 
universes showed that in an HSS universe gravitational forces alone cannot bring 
about galaxy formation. Even in the presence of an existing condensation 
(galaxy, cluster of galaxies), gravitational forces are not sufficiently powerful to 
gather large amounts of matter from regions surrounding the condensation. In 
particular, the HSS theory would require that the mass M of condensations grows 
at a rate dM/dt = 3M/T, whereas it was possible to show that only a small fraction 
of the mass of the condensation can be replenished by infall of matter during any 
period 7/3. 

This criticism cannot be directly levelled against the BGSS theory, since its 
density is not necessarily computed from the general relativistic equations and 
could be much higher than the HSS density. 

Simple calculations performed in (4) showed that gravitational forces could 
play a dominant role in BGSS galaxy formation only if the density of the universe 
was greater than a certain minimum pmin. ‘The purpose of the present paper is to 
tighten the previously used argument and thus increase the value pmin more nearly 
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to its actual value. At the same time, we will investigate just how a steady state 
galaxy may be defined. 

2. Motion of gas in the neighbourhood of a galaxy.—Consider an element of gas a 
distance r from the central condensation. At any given time, momentum con- 
servation requires that 

(p-+ Bei +780) — pf =pibt +flp = - af t+ TE ot+ 7 bp (1) 
where the mass of the galaxy, M, now contains me gravitational constant y and dp 
is the mass newly created within unit volume during time interval 5¢. The first 
term on the right is the gravitational attraction of the galaxy, the second term is the 
cosmic expansion force, and the last one represents momentum transfer from 
newly created maiter. Setting the instantaneous density p=ep, and using the 
creation rate dp = (3p,/7')dt, one can rewrite (1) as 

‘ Mr r 

r= - Stat wrt. (2) 
Here both M and « are functions of rand t. The mass of the galaxy as seen by an 
element at distance r is expressed as 


M()=M(o)+ © J.” () dt— (19 | (3) 


where M(o) is the mass at time ¢=o when the particle was at a position 7;. The 
term under the integral sign represents mass gain through continued creation of 
matter within a sphere described by the radial position of the particle; the last 
term gives the mass that would be found in the volume swept out by elements 
moving from 1; to r, if the volume were filled with gas at the general steady state 
density py). ‘This last term is necessary since we are only concerned with ‘ gravita- 
tionally active’ mass, i.e. mass in excess of the normal steady state mass pyV 
enclosed in a volume /. 

To determine « we make use of equation (3) and take the mass enclosed by two 
spherical shells at slightly different radii r, and r, 


t 
M,—M,=M,(0)~ M,(0) + { [. (nt 138) SP dt —pullrse—r)— (ri? —ra8) 


(4) 
We now choose the same initial boundary conditions used in (4), namely we assume 
that at t =o gas at a distance r; is receding from the galaxy with a speed 7;/7,, and the 
density outside the central condensation i is pp everywhere. Then M ‘(0)= M,(0). 
Letting r, —r.= Ar and r,+7r,= 2r, with Ar <r, equation (4) reduces to 





€Po = eis weet tp Po= pa I, (Ar)ridt + ar 7 mt (4a) 

To simplify the notation we define a distance ry by 
M(0)=277,°/8T?. (5) 
(This equation corresponds to equation (5) of (4).) Wealso define a constant K as 
K= 4nT%po/3. (6) 


We shall express the initial position 7; of a given element of matter by 7;=mry. Its 
position at any subsequent time ¢ will be called amr, or ar;. In this notation 


B= Ar/Ar,=1,[Aa/Ar;] + (7) 
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and equation (4a) can be rewritten as 


et 


ser) 284 ap ™ 
To integrate the equation of motion (2), we multiply both sides by r. Sub- 


stituting the expression (3) for M in equation (2) and integrating each term on the 
right side of (2) separately, we obtain, in “ new me of (5) and (6): 


+ [ale HIG) Ree SLL 


en wa 
+ 7 + — 74 (9) 
3. Mass growth.—Consider the mass enclosed by a spherical shell at distance «. 
The mass growth is AM=M(t)— M(o). From equations (3), (5) and (6), this 


eEe= 


aBdt+ — 


becomes 
AM _ 8Kn'T 3 23 
M7 a dt — (a3 |. (10) 
Suppose one requires that equation (10) satisfies the condition 
a 74>} (11) 


at time t= NT/3, where A >1 is a constant and t>7. This equation has the 
asymptotic solution 
a= aA? (12) 
where ‘a’ is constant. Substituting expression (12) into equation (10) one finds 
that (for large NV, which means «> 1) 
AM  8n'aKT 1 , 
— —1 A. I 
M27 \iogd ‘| (13) 
Through proper choice of parameters n and a, this expression can be made equal 
to any constant multiple of A’. 
One still has to satisfy oneself that the solution (12) also obeys the equation of 
motion (g), at least asymptotically. This, in fact, is the case. With 
a—-aA'T, ~4->(log A/T)«, (14) 
«> (log A)“ 





one can rewrite (9g) as 


2 
(“8 a? —> s a? — (6g 4) a So Ks LI, a at |i 5 +> 
T T? 09% 


6(log A)? f ‘ 6(log A)® P . 
+ - + t 
+ 73 P a? dt Ts : a? di (15) 
where only terms proportional to « are kept. This leads to 
(:- ig) K- 4(log A)’— 3 log A—1 - 
or 
K-> 4(log A)? + log A. (17) 


In the steady state universe one is interested in the case where A equals the 
natural number e, so that K > 5*. 


* In the HSS universe, with K=}3, one has log A=} and the asymptotic mass growth rate 
oe 3t/4t, 
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Now consider a situation in which A is chosen slightly lower than e, say 
A=e(1—x) in the limit x>o. Then 


log A = loge +log (1 —x) > 1 —x—<x?/2. (18) 
One notes that in equation (13), the terms depending on A reduce to 
eYx(1—x)%. (19) 


This means that even for small x, the mass growth will eventually fall behind 
the required value, no matter how large the coefficient in equation (13) is made 
through appropriate choice of large n and a values. 

One now defines the boundary of a galaxy in the following way. Consider 
particles that move with speeds corresponding to an A value A=(1—x)e. At 
some large distance from the central condensation there will be particles which 
move faster than this, and in particular at extreme distances A becomes identical 
with e. (This implies, of course, that a condensation must be surrounded by a 
region of density below py at some distance.) One can now define the boundary of 
the galaxy as the location of the furthest element that moves according to the 
prescription A=e(1—.x). Particles moving faster lie beyond this point and do 
not belong to the Galaxy. 

Unfortunately this definition breaks down in the interesting case of primary 
concern, namely when A =e, for then the entire Universe would fall into the 
domain of the Galaxy. 

One might note that for K >5 there is no solution of (16) for which A <e. 
At first sight one might be surprised that there even exist asymptotic solutions of 
(9) for which A >e. However, the reason for this possibility is that an expansion 
faster than exp (t/7) leads to a density less than p, in the volume enclosed by a 
shell of ‘radius’ «. This enclosed volume therefore has negative ‘ gravitationally 
active ’ mass which actually repels matter at «. 

4. Stability of the solution.—To test the stability of the solution (12) one 
examines a slightly perturbed function 


a=aA!T += at+7 (20) 


where »<,, and one writes the time derivatives as 4=a +7. Substituting 
expression (20) in the equation of motion (g) subject to «>1, one obtains an 
integral equation relating 7 to 7. We are primarily interested in the sign of 7 
relative toy. If and 7 have opposite signs the solution (12) is stable ; otherwise 
not. 

To perform the indicated integrations, one can assume that 7 and 7 are slowly 
varying functions and can be replaced by suitable average values (7) and <7), 
for simplicity in integration. This can always be accomplished by choosing a 
short enough integration interval t —t,,, where t, marks the onset of the perturba- 
tion and ¢ is the instant at which 7 is evaluated. Neglecting all terms higher than 
first order in 7 and 7, one finds after lengthy but straightforward reduction that 
terms independent of 7 again lead to equation (16), while terms linear in y and 7 
give 

[142—6+2K2—2] (4) = (n) T-{ [62 + 2(K +1) —5K27] 
+8[2*—2]} + 3[2?—2]0 (21) 
where, for brevity, we have set z = log A, 6 = 8 —« in equation (7). 
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We are interested in values of A near the natural number e. When A =e, the 
second and third terms on the right of equation (21) vanish, and the equation 
reduces to 


(9)=-4 OT”. 
The solution is stable for d=e. We rewrite (21) as 
[14542 +2271] Gj)=T-dn) Bat Bs — 3] + He? 2] 


with the aid of (17). For A <e, i.e. when o < z<1, the coefficient of {7 ) is always 
positive, and that of (7) is always negative, its absolute value being at least 25/3 
times larger than the coefficient of 6. Since | @| will usually not exceed |» | except 
possibly in some small regions, the overall effect of the 6 term will be small 
compared with the 7’~!(7) term, so that over large regions the solution (12) will 
be stable for A<e. When 7 is independent of r,, 9=o and the perturbation 
always leads to instability for A >e. 

53. Discussion.—The motion of gas particles can be computed using equation 
(9). Fig. 1 illustrates the particle trajectories for K = 5, and n values 0-85, 1-0 and 
1:18. Curve c is particularly instructive. One notes an initial decrease in 
speed at time ¢=0, when the gravitational field of the galaxy M is ‘turned on’. 
As the gas flows outward at this new, reduced speed, the momentum of locally 
created matter accelerates it away from the central condensation. However, the 
rate of mass growth determined by equation (3) eventually overpowers the repul- 
sive forces due to newly created matter, and the gas is rather rapidly decelerated, 
the direction of its motion being reversed shortly after time 7. 


2 


1-0 4 








0 Ts 2T/s T t 


Fic. 1.—Particle trajectories given in terms of the dimensionless parameter x. «* 1s plotted as a 
~~ of time for n values (a) 0°85, (b) 1:00, and (c) 1°18. The density of the universe is taken 
as K=5. 

All the trajectories of Fig. 1 were computed by means of an iteration method, 
working entirely with equation (9) witha=1anda=1att=o. In order to keep 
the labour involved within reason, it was assumed that B=« in equation (7). This 
yields a fair approximation to the actual trajectories since, at short times after t =o, 
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a is only weakly dependent on7,;; even as late as t= T/3, « = 1-261 for m= 1-0 while 
%=1°314 form=1-18. At later times, when 8 does appreciably differ from «, the 
first term on the right of equation (8) dominates the value of ¢ and in this term the 
influence of 8 is somewhat limited. Furthermore ¢« appears in only two terms 
which are roughly equal, opposite in sign, and have integrands directly propor- 
tionaltoe-!. This leads to a weak e-dependence. Still, the actual values shown 
by the curves should be taken with caution. They are intended primarily 
to indicate trends in the particle motion. 

The mass values for the three curves at times when «=o, i.e. when the gas 
just starts falling back into the galaxy, are 


n Time M Expected M 
0°85 1/3 T 1°7M, 2°7M, 
1-00 2/3 T 4 M, 74M, 
1°18 1:02 T 16 M, 21°5 M, 


One notes that the last value, i.e. 16 Mo, is an appreciable fraction of the ideally 
expected 21-5 M,, based on a simple exponential growth expectation. It appears 
that at time 7 the transient effects are becoming smaller, and one approaches the 
mass growth rates expected for K=5. At the same time, one notes that if one 
wishes the Galaxy to grow sufficiently within just a few periods 7/3, K should be 
chosen higher than 5. 

This numerical work substantiates the previously derived conclusions based 
on the asymptotic behaviour of «: the value K=5 represents a minimum density 
that would be needed if galaxies in a Bondi—Gold universe are to be formed by 
gravitational means. K-=5 denotes a density ten times that initially postulated 
in Hoyle’s relativistic steady state model. For T7=4x 10!"sec, the minimum 
density is found to be gx 10-*°gcm~*. This is considerably higher than the 
normally estimated maximum density of the universe, but the next few years 
should definitely settle on observational grounds whether or not the actual 
density of the universe is too small to sustain gravitational galaxy formation in a 
Bondi-—Gold steady state universe. 
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M.N., 122, No. 1, 1961: 
M. Harwit, Can gravitational forces alone account for galaxy formation in a 
steady-state universe ? 


p. 48, equation (5), 
for romax = (8My/3H?)!* read romax = (8My/27H’)'* 
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THE DETERMINATION OF THE INCIDENT FLUX 
OF RADIO-METEORS 


II. Sporapic METEORS 
T. R. Katser 
(Received 1961 July 26) 


Summary 


The theory relating the radio-echo rate to the incident meteor flux is 
extended to the case of sporadic meteors. It is shown that the relation 
between total influx and the echo rate observed at the equator with a suitably 
oriented aerial should be substantially independent of the distribution of 
sporadic radiants over the celestial sphere. The meteoric influx as a function 
of limiting magnitude (to +11™) has been determined from experimental 
data. Results for sporadic meteors and for several showers are presented. 





1. Introduction.—In Part I* of this work (1) formulae and graphical procedures 
were developed which enable the radio echo rate obtained with an equipment of 
given sensitivity and aerial configuration to be related to the absolute incident flux 
of shower meteors with a well defined radiant point. In extending the analysis to 
the case of sporadic meteors we must take account of the fact that the observed 
echo-rate is dependent upon the distribution of sporadic meteor radiants over the 
celestial sphere as well as upon other factors listed in the summary of Part I. 

The incident flux of sporadic meteors whose radiants lie within an element of 
the celestial sphere of solid angle dQ can be represented by (2), 


(xz) da, dQ = flux producing zenithal electron line densities between «, and 
az + daz, 
((x%,)dQ  =flux producing zenithal line densities greater than «;. 
‘These are the fluxes across a plane normal to the meteor paths. ‘The mean, over 
the whole earth, of the sporadic meteor flux through a horizontal plane is therefore 


©,(a2)=4 f, (xz) dQ (1) 


where the integral is taken over the whole celestial sphere. ‘The zenithal electron 
line density, «z, is defined as the maximum number of free electrons per unit 
height produced along the meteor trajectory ; it is related to the radio-magnitude, 
M, by (3) 


M = — 2:5 logy %z + 35 (2) 
and to the maximum electron line density, «,,, along the trail by 
tz = py (cos i (3) 


where x is the radiant zenith distance. 


* See corrections to Part I at the end of this paper. 
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Since the radio reflection points on the meteor trails are restricted to a narrow 
height range about a mean height h we can, for many purposes, consider the echoes 
as originating on the echo surface defined as the spherical surface at height h. 
From a knowledge of the equipment parameters, the minimum detectable electron 
line density, x, can be calculated for each point on this surface and a set of sensitivity 
contours plotted for a range of values of the sensitivity factor p, as in Fig. 2, Part I*. 

Some further concepts from Part I which will be valuable in the following are: 
(a) The echo plane ; this is defined, for a given radiant point, as the plane normal to 
the radiant direction and passing through the location of the observer. (6) The 
echo line, which is the intersection between the echo plane and the echo surface. 

2. Sporadic meteors observed in a conical element, dw, of the aerial beam.— 
For a given meteor distribution function, ©, we can derive the rate of radio echoes 
within a conical element of the aerial beam by arguments similar to those used in a 
previous work concerned with the meteor height distribution obtained from radio- 
echo observations (4). The geometry is represented in Fig. 1. The conical 
element has half angle 4 and the radiants of the observed meteors consequently 
lie in the shaded strip of the celestial sphere of width 2. The elevation of the 


Zenith 











Fic. 1.—Meteors observed in the conical element of the aerial beam (with semi-angle ) have 
radiant points Q (co-ordinates (€, £)) within the shaded strip of the celestial sphere. 


element above the horizontal is ¢, measured at the location of the observer, and ¢ 
where it intersects the echo surface. 

Let us consider meteors radiating from the point Q with co-ordinates (£, ¢) 
(see Fig. 1). The echo plane corresponding to Q will intersect the conical beam 
element, the angular intercept being 2(? — ¢?)"* and hence, referring to the echo 
plane geometry depicted in Fig. 1 of Part I, the linear intercept between the conical 
element and the echo line is 2R(s? — €?)"* (cos 0)-! where R is the range to the echo 
surface and @ is the angle between the beam element and the normal to the echo 


*p is inversely proportional to the minimum detectable line density and may be normalized to 
a maximum value of unity. 
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line (in the echo plane). We can now use equation (14) of Part I, which relates 


‘OW . . . . 
the rate of echoes observed over unit length of the echo line to the incident flux 
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7 of meteors, and obtain the echo rate observed within the beam element due to 
ron meteors radiating from an element of solid angle dé d¢ on the celestial sphere : 
yity 2_ 72)/2 
I*, Rate from df af = See ( . ) dé dé (4) 
ie sin x cos 0 cos x 
lto where H = atmospheric scale height, 
lhe Tis defined by equation (15) of Part I, 
r x = zenith distance of Q, measured on the echo line, 
aa a = minimum detectable line density. 
vies Integrating (4) with respect to ¢ between the limits —% to + and putting 
ef sin y cos#=sin¢ we obtain: 
ica 
tly Rate from strip €—(€+dé) = a o( = ) dé. (5) 
the “sing — \cosx 
Putting ms? = dw and integrating with respect to €, we get: 
[2 
Rate within beam element dw = smal! 0) ( _ ) ae ] dw. (6) 
singL J 2/2 \cosx 
3. Total rate per unit range interval.—Consider a beam element of 
rectangular cross section at elevation ¢) and azimuth # and hence of solid angle 
dw =cos $y d¢,d%. The range interval on the echo surface intercepted between 
elevations ¢) and (¢)+ dd) is dR= R(tan¢)—1 d¢dy, hence 
Rdw _ 8b apag. 
. sind cos¢ 
Now from equation (32) of Part I 
cosdy h 
cos¢ Ry 
where R, is the earth’s radius. To a sufficient approximation h/R, can be 
neglected, compared with unity, and hence from (6) we get: 
Total rate between ranges R and (R+4dR) and azimuths # and (# + d#) 
7/2 “ 
= HIARas | o(=) dé. (7) 
—n/2 Cos x 
id The distribution function @ may be represented by 
O(az)oc az), 
¢) Now cos y=cos¢cosé~cos¢,cosé and p=%/%; % is the minimum detectable 
“a line density at range R in the direction of the beam axis and ¢y is the elevation to 
“ range R on the echo surface. Using (7) we thus obtain the total echo rate NdR, 
“ observed between ranges R and R+dR: 
) 


NdR=HI ar {{o ( 


=) (pcos) 1dedo (8) 

















268 T. R. Kaiser, The determination of Vol. 123 


where the integral with respect to # is taken across the beam at constant range R. 
If the aerial beam is sufficiently narrow in azimuth, so that © is reasonably constant 
over a distance on the celestial sphere of the order of the beam width, (8) becomes 


NdR~ H1AR | p s— rao [™ o( a ) eos dé. 
—2/2 COS Pp ode (9) 


The integral with respect to # can be evaluated graphically, using the sensitivity 
contours projected onto the echo surface (Fig. 2 of Part 1). The difficulty lies 
in the second integral and the absolute incident flux cannot in general be deduced 
directly from the echo rate obtained at a single station, unless the form of the 
radiant distribution and the exponent sare known. ‘The latter can be determined 
in several ways (Browne et al. (5), Kaiser (2, 4, 6)). Values previously obtained 
are near to 2:0, while recent determinations at Sheffield (7, 8) suggest a small 
diurnal variation in s, with a mean value of 2-17 for meteors between 8th and 11th 
magnitude. 

4. The radio-echo rate at the equator.—'Vhe case of an aerial beam, narrow in 
azimuth, and directed eastwards or westwards at the equator, is important in that 
the radiants of the observed meteors lie along a great circle from the north to the 
south pole and in the course of a day this great circle sweeps the whole celestial 
sphere. Let 7 radians be the hour angle of this great circle and let N dR be the 
mean diurnal echo rate observed between ranges Rand R+dR. Then, from (9) 
we obtain, 


2x NdR= HIAR | p - do | nf o(= ) eos" lé dé dT. (10) 
—n/2 cos do 


Now dQ=cos€dé dT is an element of solid angle on the celestial sphere, hence 


an NdR=HIAR | p*-140 | o( to ) cost-*8 d0. (11) 
COS po 
Since there are relatively few sporadic meteors with radiants near the poles 

(i.e. with € near go“), negligible error will be involved if we put cos*-*é = 1, hence 














Equation (12) can be extended to the case where the echo rate is evaluated over 
a finite range interyal by determining fp*-!d#@ as a function of R (using the 
sensitivity contours of Fig. 2 of Part 1) and integrating graphically. 

The importance of (12) is that it enables the absolute mean influx, 0,, of 
meteors over a sphere at the Earth’s distance from the Sun to be deduced from 
the mean diurnal rate, N, irrespective of the form of the distribution of sporadic 
meteor radiants over the celestial sphere. Thus, for instance, the seasonal 
variation in the influx of meteors brighter than a given zenithal magnitude can be 
determined. In addition the results are directly comparable with satellite 
observations which determine the influx of meteoroids over a sphere as a function 
of their momenta. 

5. Some preliminary results.—Since equatorial measurements of sporadic 
radio-meteor rates have not been made, we are obliged to use data obtained at 
higher latitudes. A recent analysis (7, 8) has been made using sporadic meteor 
rates obtained with a 17 Mc/s radio echo equipment at Sheffield (latitude 55°N). 
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The aerial was directed 65°W of N, at low elevation (sensitivity contours are 
given in Fig. 2 of Part I), and it is considered that the total daily rate during the 
autumnal equinox (when the apex of the Earth’s way is at maximum declination) 
should approach the equatorial value. ‘The observed mean rate between ranges 
300 and 600km, during September 1958, was 500 per hour. ‘The limiting radio 
magnitude was 10:8 +0°5 


(=.= Wo =5x ro"em ) 


and the mean influx, ©,, was evaluated, using (12), as 1-9 x 10-!”m~*sec~! which 
corresponds to 8-2 x 10° over the whole Earth per day. 

The theory of Part I has been used to deduce the magnitude distributions and 
incident fluxes in the Geminid and Arietid showers using observations made at 
Sheffield. The results are illustrated in Fig. 2 which gives the shower and sporadic 
fluxes as a function of limiting magnitude. Some previously reported data 
(2, 5, 6) are also included. ‘To facilitate direct comparison between the shower 
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Fic. 2.—@, m=? sec-1=mean incident flux over the whole Earth of meteors with zenithal 
electron line densities greater than az (brighter than magnitude M). 


Sporadic meteors (Sheffield, 1958 (7, 8)) 
Arietids (Sheffield, 1959) 

Arietids (Jodrell Bank data (2, 6)) 
Geminids (Sheffield, 1958) 

Geminids (Jodrell Bank data (2, 6)) 
Perseids (Jodrell Bank data (2, 6)). 
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and sporadic fluxes, the values for the former are given as the mean influx over a 
sphere (which is one quarter the flux through a plane normal to the radiant 
direction). 

6. Discussion.—An important conclusion from the present work is that the 
relation between the total sporadic meteor influx and the radio echo rate observed 
with a suitably oriented aerial on or near the Earth’s equator is substantially 
independent of the distribution of sporadic meteor radiants. ‘This makes it 
possible to determine the variation of activity with solar longitude by making 
regular equatorial observations over a period of a year. ‘The value of such an 
observational programme is further enhanced in view of the accumulation of 
related data concerning the influx of impinging meteoroids on earth satellites. 

An interesting feature of Fig. 2 is that the faint shower meteors are swamped 
by the sporadic background. ‘This is true of the Arietids over the whole range 
of magnitudes covered by the observations, while the Perseid and Geminid 
fluxes exceed the sporadic level only for meteors brighter than about magnitude 
+5, (a similar result has been reported by Kresak and Vozarova (9g) from visual 
observations). ‘These shower fluxes refer only to the period of maximum shower 
activity, whereas it is known that in general the magnitude distribution in a 
particular shower changes during the Earth’s passage through the meteor stream 
(5). ‘These effects are particularly important in connection with problems of the 
evolution of meteor streams and further data are needed. ‘The theory of Part I 
should prove of value for their quantitative reduction. 

Some of the limitations of the theory have been discussed in the previous 
paper (1). One simplifying assumption is that the maximum meteor echo ampli- 
tude depends only upon the electron line density in the vicinity of the radio 
reflection point and on the equipment parameters. It has been shown (6) that 
rapid diffusion may cause a reduction in amplitude, the effect being greater the 
greater the height and the shorter the radio wavelength; a similar reduction is to 
be expected due to the finite initial radii of the trails. "The consequent reduction 
to be expected in the echo rate has been discussed by Peregudov (10) and Fialko 
(11, 12), while Greenhow and Hall (13) have recently presented experimental 
evidence for a significant initial radius effect. ‘These limitations may be minimized 
by judicious choice of equipment parameters, particularly the wavelength, but 
their quantitative effect on the observed echo rate certainly warrants further 
careful investigation. 


Acknowledgments.—I am indebted to my colleagues J. P. Baldwin and D. 
Brooks for their co-operation in obtaining and reducing the observations presented 
in Section 5, and to L. M. G. Poole for pointing out some errors in the previous 
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Corrections to Part I (M.N., 121, 284, 1960). 


Equations (1), (2), (3) and (15) in Part I should read 
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The last equation also needs correction where it appears in the title of Fig. 7. 
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Summary 


Radial velocities and UBV colours are presented for 21 newly identified 
southern subdwarfs. A redetermination of the solar motion relative to the 
subdwarfs using the new material gives a velocity 131+27 km/sec in the 
direction of galactic rotation. There is little difference between the kinematic 
properties of subdwarfs with high, and those with low ultraviolet excesses. 





Introduction.—Spectroscopic subdwarfs have mostly been found as a by- 
product of routine programmes of radial velocity observations of proper motion 
stars. In view of the scarcity of observatories in the southern hemisphere, and the 
pressure there of more immediately interesting problems, such work has normally 
been carried out from northern observatories only. This is apparent in Eggen 
and Sandage’s useful list (1) of known spectroscopically identified subdwarfs, the 
distribution of which with declination is rather strongly non-uniform. 

The subdwarfs, however, are mong the most easily identifiable representatives 
of the ‘‘ Extreme Population I1’’ component of the Galaxy in the general field, 
and it is desirable that as many of them as possible should be found. During a 
recent visit to the Radcliffe Observatory, Pretoria, South Africa, the author con- 
ducted a spectroscopic search for southern subdwarfs among likely stars of high 
proper motion. As a finding list for these, Luyten’s useful catalogue of southern 
proper motion stars (2) was used, referred to as L'I'T in what follows. While this 
is not the most desirable method of selecting low luminosity stars, since it biasses 
the selection in favour of stars with high velocities relative to the Sun, it is the only 
one available, as there are, at present, no suitable lists from work based on objective 
prism or similar surveys which would not suffer from this selection effect. 

Observations.— Spectrograms of selected stars were taken with the ‘‘d’’ camera 
of the two prism Cassegrain spectrograph attached to the 74-inch Radcliffe 
reflector (dispersion 86 A/mm at Hy). Normally, a first exploratory spectrogram 
was taken at a slit width of o-20mm, which gives a projected slit width on the 
plate of 27 microns, equivalent to about 160 km/sec at Hy. Further plates of 
stars with interesting spectra or a high radial velocity were taken at a slit width 
of o-14mm so as to produce purer spectra for radial velocity measurements. 
Kodak 103aO emulsion was used throughout. 

An attempt was made to classify all the spectra from plates of standard stars in 
the Radcliffe Observatory’s plate library. Joy (3) has given a definitive description 
of the spectral peculiarities shown by subdwarfs. The most general feature is a 
marked weakening of all the metallic absorption lines. Joy’s description has to be 
interpreted in the light of the current belief that most of the stars which were 
formerly classified as A-type subdwarfs ought properly to be classified F, following 
analyses of high resolution spectra of some of these (Chamberlain and 
Aller (4), Baschek (5), Aller and Greenstein (6)). Furthermore, there is 
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some doubt as to the existence of spectral differences between K or M dwarfs 
and subdwarfs; Unsold (7) believes on theoretical grounds that cool subdwarfs 
should show no spectral peculiarities and disagrees with Joy’s reported observa- 
tions of metal weakening in some M-type stars. 

‘The present results are confined to F- and G-type stars, where the spectro- 
scopic effects are certain, and here weakness of the metallic absorption lines relative 
to the hydrogen lines is the sole criterion for the selection of subdwarfs. The degree 
of metal weakening necessary before a star is called a subdwarf is a rather arbitrary 
quantity. Here the term is restricted to stars showing a quite large deviation from 
normality, where the use of an MK type, even if moditied by a ‘‘p’’, no longer 
gives even an approximate description of the spectrum. In these stars, all the 
metallic lines, and the G band of CH are very much weakened relative to the 
hydrogen lines which, together with the H and K lines of Catt are frequently 
the only lines visible. It is most undesirable to apply an MK classification to a 
star which is recognisably far from normal, since this can lead to considerable 
confusion in the interpretation of the results. This difficulty is partly avoided, 
but not entirely overcome, in the present list by the deliberate use of unsubdivided 
spectral classes to describe the stars. ‘Thus a type ‘‘sdF’’ implies a general 
resemblance to an F-type spectrum, although it may be applied to a star which on 
other evidence seems to have the effective temperature of a normal late-type G 
star. Beyond refusing to give a spectral class at all (which is possibly the only 
final answer to the problem), this procedure seems the best that can be followed at 
present. 

Each spectrogram was measured for radial velocity using the Casella projection 
measuring machine at the Radcliffe Observatory. A set of iron arc wavelengths 
was used which was substantially the same as one found by Kinman (8) to be 
suitable for use with this camera at this dispersion. All stellar lines in the spectra 
of subdwarfs are weak, and there are frequently only a few lines suitable for radial 
velocity measurements: Hy, Ca 1 4226, Hd, Fei 4045, 4063, 4071, Catt K are the 
most easily measurable. Since the wavelength of an absorption feature is the 
weighted mean of the wavelengths of its contributors, and since the relative 
weighting will be different in stars of low metal content, Petrie’s wavelengths (9) 
derived for normal stars are not in general applicable to the subdwarfs. The 
present velocities were reduced using laboratory wavelengths for all the lines. An 
examination of the residuals from each line showed that no corrections to individual 
lines were necessary, and none were applied. Inany case, the systematic difference 
in radial velocity between two possible systems occurring through the use of 
different wavelengths is much smaller in the present case than the average observa- 
tional error of each determination, since only a comparatively few lines and plates 
have been available for measurement. 

In forming the mean velocity for each star from individual plates, determina- 
tions with standard errors (from the internal agreement of velocities from different 
lines) up to 8-5km/sec were given unit weight, determinations with standard 
errors between 8-5 and 15 km/sec were given half weight, and determinations with 
standard errors greater than this were rejected. These figures were suggested 
from an examination of the distribution of the individual standard errors from 
each plate. ‘The average standard error of a final radial velocity is 3-5 km/sec. 

During his visits to South Africa in 1960, Dr O. J. Eggen of the Royal 
Greenwich Observatory, Herstmonceux, measured U, B, and V for most of the 





new 
or t 
alloy 


follc 
have 
star: 
star: 
belc 
wot 
dw 
mos 
aris 
SCO] 
see 
Pog 
pro 
dw: 
clus 


ma; 
and 


“ce 


$l 


are 
onl 
i.€. 
is | 
evi 
altl 
Iti 
ma 
glo 


nu 
val 
sys 
F 


the 
en: 
ex’ 


10 
in 
of 
dif 


1. 123 


1warfs 
lwarfs 
Servas 


ectro- 
‘lative 
legree 
itrary 
| from 
onger 
Il the 
to the 
1ently 
n toa 
srable 
vided, 
vided 
-neral 
ch on 
pe G 

only 
ved at 


ction 
ngths 
to be 
ectra 
-adial 
‘e the 
s the 
lative 
s (9) 
The 

An 
idual 
‘ence 
se of 
2rva- 
lates 


1ina- 
srent 
dard 
with 
sted 
from 
/sec. 
oyal 
F the 





XUM 


No. 3, 1961 Some new southern subdwarfs 275 


new subdwarfs with the Radcliffe 74-inch reflector and with the 26-inch refractor 
or the 18-inch reflector at the Cape. I am greatly indebted to Dr Eggen for 
allowing me to use his results. 

Are these stars true subdwarfs ? Sandage and Eggen (10) and Melbourne (11), 
following a suggestion originally due to Schwarzschild, Searle and Howard (12), 
have shown that the relation between effective temperature and B—V for dwarf 
stars depends upon the composition of the star, and they suggest that some 
stars which appear as subdwarfs in, say, the (/,, B—V) diagram, would not lie 
below the normal main sequence in an (V/,,,, log T,) diagram, and in such a diagram 
would not be called subdwarfs. ‘The usual method of identification of ‘‘sub- 
dwarfs ’’ from their spectral peculiarities picks out just those stars which deviate 
most from the normal in their (B—V, log 7) relation. ‘The question therefore 
arises of whether they should still be called subdwarfs. Since these spectro- 
scopically peculiar stars, more than any other selection of possible ‘‘ subdwarfs ”’ 
seem most closely connected with globular clusters and other representatives of 
Population II, which have similar kinematic, spectroscopic and photometric 
properties, it would seem that a better name for them would be ‘‘ Population II 
dwarfs’’. This idea is strengthened by the fact that the main sequences of globular 
clusters in the (M,, B—V) diagram lie one or two magnitudes below the normal 
main sequence when the RR Lyrae stars in the clusters are fitted to an absolute 
magnitude of about o™-5. However, in order to avoid unnecessary confusion, 
and since these stars probably are subdwarfs in some HR diagrams, the usual term 
“‘subdwarf’’ has been retained here. 

Finally, one should emphasise that normal spectroscopic luminosity criteria 
are not present in these stars, since all the metal lines are so weak. Almost the 
only information we have from the spectra concerns the degree of metal weakness, 
i.e. (perhaps) the ‘‘ population’’ type. The supposition that the stars are dwarfs 
is based on their high proper motions alone. In particular there is no direct 
evidence to indicate that these stars lie below the ‘‘normal’’ main sequence 
although indirect evidence (e.g. the ultraviolet excesses) suggests that they do. 
It is just possible that some of the stars could be intrinsically brighter than normal 
main sequence stars (perhaps like those stars just above the ‘‘turn-off’’ point in 
globular clusters) although, in view of the high proper motions, this is unlikely. 

Results.—Table I gives the results of radial velocity and photometric measure- 
ments for the new southern subdwarfs, which are identified by their LTT 
number. ‘The values of V, B—V, and U— Bare by Dr O. J. Eggen. Where the 
value of U— B is given in brackets, it is a measure on the Cape refractor (U — B), 
system. ‘These new stars provide a useful complement to the previously known 
F and G subdwarfs in Eggen and Sandage’s list (1), and it has therefore seemed 
worthwhile to use the new data in a redetermination of the solar motion relative to 
the subdwarfs. The restriction to F- and G-type spectroscopic subdwarfs 
ensures homogeneity of the data, and provides probably as pure a sample of 
extreme Population II dwarfs as is available at present. 

Eggen and Sandage’s subdwarfs, plus the new southern ones, give a total of 
103 stars for analysis. Preliminary least-squares solutions for the solar motion 
included a K term in the analysis to represent a possible expansion or contraction 
of the system of subdwarfs. However this was found to be not significantly 
different from zero, and the solution was repeated omitting the K term. For 














276 T. F. Deeming Vol. 123 


TABLE I 

Radial velocities, spectral types and photoelectric colours for 21 new subdwarfs 
LTT Spectrum V B-V U-B R.V (km/sec) n 

877 sdF-G +227+8 (5) 
2237 sdF 10°63 +0°54 (+ 1°50) +320+6 (5) 
2679 sdF 9°05 +0°46 (+1°56) + 71+2 (3) 
2744  ~sdG +247+5 (3) 
3483 sdF-G IO"IO +0°48 — 0°07 —110+5 (4) 
3774 sdF I0°OI +0°43 —0°!3 +194+2 (5) 
4226 sdF 8-96: +0°47 —o'2I +2224+3 (5) 
4519 sdF-G 10°39 +0°44 +121+2 (3) 
4859 sdF 9°87 +0°47 —o'16 +220+5 (5) 
5145 sdF + 14443 (3) 
5468 sdF 8-68 +0°48 —o'18 — 23+6 (5) 
7043 sdF + 10443 (4) 
7182 sdF-G 10°22 +0°57 —o'19g +364+2 (6) 
7321 sdF 10°92 +0'50 —Or'ls +1064+3 (5) 
7662 sdF 10°72 +0°50 —o0'28 ia 
7668 sdF 9°62 +0°44 —o0'28 + 1727 (2) 
7862 sdF 10°12 +0°42 —o'26 —146+4 (3) 
8578 sdF 10°52 +0°40 (1°45) +112+1 (5) 
9074 sdF 9°64 +040 — 0°24 + 47+3 (4) 
g201 sdF 10°05 +0o°41 —0°25 — 83+7 (5) 
9509 sdF IO'I5 +0°66 —0'05 + 72+4 (4) 


103 stars, the mean subdwarf motion relative to the Sun, in equatorial components 
(epoch 1950), was found to be 

X= —67+25km/sec 

Y= +84+27km/sec 

Z= — 103 + 27km/sec. 
(For a discussion of the application of the least-squares method to this problem, 
see ‘T'rumpler and Weaver (13), p. 181). On transforming to galactic velocity 
components U, V, W, where U is directly away from the galactic centre, W in the 
direction of the north galactic pole, and V perpendicular to U and Win the direction 
of galactic rotation, and using the solar motion determined by Woolley (14), we 
find the mean subdwarf velocity relative to the local standard of rest to be: 

U= -—g9+27km/sec 

' V = —131+27km/sec 

W= +3+26km/sec. 

There is thus no significant deviation of the vertex of the subdwarf motion 
from the direction of galactic rotation. Ina similar analysis of the radial velocities 
of 79 subdwarfs, Fricke (15) found a velocity of the Sun relative to the subdwarfs 
of 148+ 28km/sec directed toward galactic longitude /'=75° + 8° (old system). 
Applying Woolley’s value for the solar motion relative to the local standard of rest, 
this becomes 135 + 28 km/sec directed toward 1=77°. ‘The agreement between 
the two velocities is good, but Fricke’s result shows a deviation of 19° in the vertex 
of the motion from the direction of galactic rotation. This may be due to his 
inclusion of some trigonometrical subdwarfs based on rather uncertain parallaxes, 
as well as some high velocity stars, whose subdwarf nature must be in doubt. 
The present results indicate a mean velocity dispersion of 135 km/sec; Fricke’s 
value was 110 km/sec. 
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Further information about the subdwarf motions may be obtained from the 
proper motions and parallaxes. For most of the known F and G subdwarfs, 
trigonometrical parallaxes are not available. However the blanketing theory of 
subdwarfs (10) gives a method of estimating the absolute magnitude of a sub- 
dwarf from its observed B— V colour, and its ultraviolet excess, 6(U—B). Hence 
we can derive a ‘‘ photometric parallax’’. ‘The accuracy of these parallaxes 
depends on the accuracy of the measurements of U, B, and V, and on the correct- 
ness of the table of blanketing corrections given by Sandage and Eggen (10). Of 
these two, the latter is almost certainly much the greater source of error. There is 
not yet sufficient information to give a judgement on the correctness of this, but 
we may tentatively suppose that an error not greater than +0™-5 in absolute 
magnitude may be expected, giving a relative error in parallax, 57/7, of 0-23, anda 
similar error in the derived transverse velocity component. 

The proper motions (mostly from the extensive publications of the Cape 
Observatory), the radial velocities, and the photometric parallaxes were 
combined to compute space velocities in galactic velocity components for each star 
which has been measured on the UBV system, and which has an ultra-violet 
excess greater thano™-10. The stars were then divided into two groups according 
to ultraviolet excess: (a) those with ultraviolet excess between o™-10 and o™-21, 
and (4) those with ultraviolet excess greater than o™-21, and mean motions and 
dispersions in each velocity component determined for each group. ‘The results 
were: 


(a) Low excess stars 


U=— 29+22km/sec of, =105+15km/sec 

V=-—113+20km/sec oy, = 93+14km/sec 

W=+ 12+12km/sec og = 55+ 8km/sec 
(6) High excess stars 

U=— 11+22km/sec Oy =100+ 16 km/sec 


HT 


V=—141+17km/sec oy 
W=+ 2+18km/sec op 


75 +12km/sec 
78 + 12 km/sec. 


Note that the velocity dispersions, o,, oy, ow are not corrected for possible 
errors in the adopted photometric parallaxes, since these are probably systematic 
rather than random errors. There is apparently a slight deviation of the vertex 
of subdwarf motion for the stars with low ultraviolet excess which is in the same 
direction as that found by Fricke in his analysis of subdwarf motions, but it is 
hardly significant. The high excess subdwarfs show no evidence of a deviation of 
the vertex of the motion. There is only a slight difference of 28 + 26 km/sec in 
mean motion between the low excess stars and the high excess stars. There is a 
significant difference in the velocity dispersion in the W direction between the low 
excess and the high excess stars, and this implies a different distribution of stars 
with height above the galactic plane. The low excess stars with the smaller 
dispersion of W velocities must extend less far above the galactic plane than the 
stars with the higher velocity dispersion. 

It is interesting to compare these results with similar results derived by Kinman 
(16) for the globular clusters and RR Lyrae stars. He finds a mean solar motion 
relative to the globular clusters of 167 + 30km/sec, similar to that for RR Lyrae 
stars with periods between 04-5 and 04-6. This is just in agreement, within the 
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limits of the observational errors, with the value for all the subdwarfs. The mean 
velocity of all the ‘‘ Population II’’ objects—g\obular clusters, RR Lyrae stars, 
and subdwarfs—is then 131 + 18km/sec, but this value has no particular astro- 
physical significance unless one has grounds for supposing that all the Population II 
objects had an almost identical galactic origin, in which case their kinematic 
properties should be the same. Some of Eggen’s moving groups of subdwarfs (1) 
suggest that this may be the case, but the evidence is by no means conclusive. 

Note that if the commonly accepted value for the circular velocity of Population 
I objects at the Sun, of about 220 km/sec, is true, then the Population II stars have, 
on average, a quite appreciable galactic ‘‘ rotation ’’ of about go km/sec at the Sun. 
The circular velocity is only very imperfectly known, and it might be possible to 
accept a smaller value than 220km/sec. In this case, however, we should find 
that several of the subdwarfs have retrograde orbits, and it can be argued that a 
retrograde orbit is a priori improbable, since a very close (and statistically 
improbable) stellar collision is required to reverse the direction of a star’s orbit, 
and a star is unlikely to have been formed initially with a retrograde galactic orbit, 
unless there was a very high degree of turbulence in the protogalactic gas cloud, of 
the order of 100 km/sec, r.m.s. velocity. It will be interesting to await the develop- 
ment of proper motion systems based on the extragalactic nebulae, since this seems 
to offer the most promising line of approach for a precise determination of the 
local circular velocity. 
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THE RADIO EMISSION FROM NORMAL GALAXIES 


III]. OBSERVATIONS OF IRREGULAR AND EARLY-TYPE GALAXIES AT 158 MC/S 
AND A GENERAL DISCUSSION OF THE RESULTS 


R. Hanbury Brown and C. Hazard 
(Received 1961 August 31)* 


Summary 


A survey of the radio emission at 158 Mc/s has been made from twenty 
bright spiral, four irregular and four early-type galaxies. ‘This paper reports 
the measurements of the irregular and early-type galaxies; the results 
obtained from the spirals were given in the two previous parts. ‘The variation 
of the radio index R with type of galaxy is discussed; R is defined by the 
equation R=m,—my, where m, and my, are the radio and photographic 
magnitude respectively. It is shown that the value of R for irregular galaxies 
is significantly greater than that for Sb and Sc spirals, and that this is probably 
also true of early-type galaxies. 





1. Introduction 

The present paper completes an account of a survey of bright normal galaxies 
which was carried out in 1958 with the 250 ft steerable paraboloid at Jodrell Bank 
working at 158 Mc/s. <A normal galaxy, following Baade and Minkowski (1), is 
taken to be a normal member of the population of galaxies observed photographic- 
ally. In Paper I (2) we described the apparatus and also the results obtained on 
M31 and M33; in Paper II (3) we reported measurements on twenty spirals of 
type Sband Sc. Inthe present paper (III) we report some observations of a few 
irregular and early-type galaxies; we also discuss how the radio emission from a 
normal galaxy depends upon type. 


2. Observations 

An attempt was made to observe the radio emission at 158 Mc's from four 
irregular galaxies of magellanic type and four early-type galaxies. A radio source 
was observed close to the position of 1C 1613; the regions of sky near the other 
galaxies were complicated by the proximity of other sources and irregularities in 
the general background radiation and it was possible only to set an upper limit to 
their emission. ‘The results are shown in Table I. We have tabulated Mg 
the photographic magnitude according to Holmberg (4) or de Vaucouleurs (5), m,, 
the observed integrated radio magnitude at 158 Mc/s, and R the radio 
index, where R=m,—m,,. We have not corrected the values of R, as discussed 
in Paper II, for the absorption of light ; an estimate of these corrections shows that 
they are less than 0-5 mag. and for the purposes of the present paper they can be 
neglected. 


* Received in original form 1961 July 10. 
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TaBLe I 


Observations at 158 Mc/s of irregular and early-type galaxies 


Galaxy Type Mog m, R=m,—Mgg 
NGC 
4449 Ir1 9°9 >10°6 > +0°7 
4214 Irt I'l SII‘o 2 +0°9 
6822 Ir 9°2 >10°7 >+1°5 
IC 1613 I 8-9* grl +02 
4594 Sa 9°2 > 10°6 >+1°4 
3368 Sa 10°! >10°4 > +03 
4472 E 9°3 > 10°2 2+0'9 
1023 SO 10°5 > 98 =—-0'7 


* The optical data for this galaxy has been taken from de Vaucouleurs (5). 


3. The radio index of different types of galaxy 


3.1. The radio index of Sb and Sc spirals.—It was shown in Paper II that the 
average radio index of Sb and Sc spirals is + 1-4 and + 1-2 respectively, the asso- 
ciated r.m.s. deviations about these mean values being +0°8 and +0-7. These 
conclusions were based on measurements of sixteen galaxies selected from a survey 
of twenty bright spirals. 

3.2. The radio index of early-type galaxies.—In view of the scarcity of observa- 
tional data we shall assume that types E, SO and Sa all have the same average radio 
index and consider these three types as a single group. 

Table II is a list of all the early-type galaxies which have been surveyed 
individually for radio emission. For the purpose of comparison the measure- 
ments at different frequencies have been scaled to 158 Mc/s on the assumption that 
the spectral index is —o0-6 since this was the median value found by Harris arid 
Roberts (10) for all sources in their survey. If also, following Harris and Roberts, 
we assume that the dispersion about this mean spectral index is such that 50 per 
cent of the sources have values in the range — 0-4 to — 0-8, then we can calculate 
the associated uncertainty in scaling the values in Table I1; the greatest 
uncertainty, which arises in scaling the 1400 Mc/s results, is about + 0-5 mag. 
and is not significant in the present discussion. 

The table shows that, out of a total of eighteen galaxies, in only one case 
(NGC 4526) was a radio source observed; the remainder were below the limit of 
detection and may therefore have any value of radio index greater than that shown 
in Table II. The results cannot, therefore, be used to derive a mean radio index, 
but they can be used to estimate a lower limit if we assume a value for the dispersion 
about the mean. ‘Thus if we assume that individual early-type galaxies have radio 
indices which are distributed about the mean with a similar dispersion to that 
observed for Sb and Sc spirals ( + 0-7 mag.), then the probability of occurrence 
of the distribution of limits given in Table II can be calculated for different values 
ofthe mean. This calculation shows that there is roughly a 50 per cent probability 
that the observed distribution would arise by chance if the mean index is as low as 
+ 2°5, and only a 10 per cent probability, if the mean value is +2-0. Thus, given 
the assumed dispersion, it is probable that the mean radio index of early-type 
galaxies is greater than + 2:0 and is therefore greater than that (+ 1-3) of Sb and 
Sc spirals. 
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TaBLeE II 
Observations of single early-type galaxies 
Galaxy Type Myg m,(158 Mc/s) R (158 Mc/s) 
NGC 
185 E 10°3 >11°0(H) > +07 
205 E 8-9 > 102 (L) > +13 
1023 so 10°5 > 98 (BH) >-0o'7 
1052 E II‘7 > 11°0 (H) >-o7 
1291* SBO 9°4 >11°1 (M) >+1°7 
3r15* E 10°2 >11°3 (M) >+1'1 
3368 Sa Io'l > 10°4 (BH) > +03 
3379* E 10°3 > 10°8 (Le) > +0°5 
3623 Sa 10°2 > 10°8 (Le) > +06 
4278 E I1I‘2 >11°0(H) >—o'2 
4382 sO 10°! > 10°0 (Le) >-o'! 
4472 E 9°3 > 10°2 (BH) >+0°9 
4494* E 10°8 > 10'1 (Le) >-o7 
the 4526* Sa 10°5 100 (Le) —0'5 
4594 Sa 9°2 >11°2 (M) > +2°0 
SSO- 4636* E 10° > 10°3 (Le) > +03 
hese 4649* E 9°9 > 10°3 (Le) > +03 
rvey 5866* SO II‘! > 11'0 (H) >-o'! 
rva- H D. S. Heeschen (6) at 1400 Mc/s. 
site L M. I. Large (7) at 408 Mc/s. 
BH Present survey at 158 Mc/s. 
M B. Y. Mills (8) at 85 Mc/s. 
yed Le P.R.R. Leslie (9) at 178 Mc/s. 
ire- * Optical data according to de Vaucouleurs (5); the data for the other galaxies 
that was taken from Holmberg (4). 
aiid : 
rts, 3.3. The radio index of irregular galaxies.—Table 111 is a list of the available 
per observations of irregular galaxies. ‘The measurements have been adjusted to 
late 158 Mc/s on the assumption that the spectral index of the sources is —0-6. All 
test the galaxies in the list contain a substantial proportion of Population I; in the 
lag. classification of Holmberg they appear as type Ir1, and in the system of 
de Vaucouleurs they are referred to as Magellanic. One galaxy, IC 1613, appears 
sase to have a radio emission which is markedly different from the others. We have 
t of therefore surveyed this galaxy also at 240 Mc/s and have confirmed the presence 
ywn 
lex, Tase III 
oe Observations of single magellanic irregular galaxies. 
“a9 Galaxy Mog m, (158 Mc/s) R (158 Mc/s) 
a LMC 06 3°3 (M) +2°7 
nee § SMC 2°5 5°4 (M) +2°9 
ues ff NGC 55 7°8 10°8 (M) +3°0 
lity f NGC 4449 9°9 > 10°6 (BH) >+0°7 
y as NGC 4214 10°! >11'0 (BH) >+0'9 
en E NGC 6822 9°2 > 10°7 (BH) >+1°5 
IC 1613 8-9 g*1 (BH) +02 
ype 
and M B. Y. Mills (8) at 85 Mc/s. 
BH The present survey at 158 Mc/s. 
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an source close to the optical position; however, the value of radio index mea- 
sured at 240 Mc/s was +1-°3 which differs considerably from the value of + 0-2 
obtained at 158 Mc/s. This discrepancy suggests that the observations are 
affected by confusion with other sources and that the association of this radio 
source with IC 1613 must be regarded as doubtful, although the possibility cannot 
be excluded that IC 1613 is an excited or radio galaxy. 

If we exclude the result on IC 1613, then the data in Table III indicate that 
the average radio index of magellanic galaxies is certainly greater than that of 
spirals; the measurements on NGC 55 and the Magellanic Clouds suggest that 
it has a value of about + 3-0 

3.4. Summary.—Fig. 1 summarizes the results of the previous sections. It 
shows that the ratio of total radio emission to the total light, as measured by the 
radio index is greatest for Sb and Sc spirals and decreases significantly for 
irregular galaxies. It also seems likely that it decreases for early-type galaxies, 
although this conclusion cannot be established by the data. These results 
confirm the conclusions reached by Mills (11) from earlier surveys. 

There appears to be no simple correlation between the radio index of normal 
galaxies and their mass, linear size or colour index. However, the results in 
Fig. 1 suggest one particular conclusion. It appears that the ratio of coronal to 
disk radiation in an irregular galaxy is much less than that found in an Sb or Sc 
spiral. ‘The absence of a corona around irregular galaxies is suggested by an 
examination of the radio isophotes of the Magellanic Cloud (8), but their large 
angular size and proximity to our own Galaxy makes this conclusion somewhat 
uncertain. However, if we assume that NGC §5 and the other four irregulars in 
Table I are similar to the Magellanic Clouds, then all the measurements taken 
together show that irregulars do not have extensive and powerful radio coronae 
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Fic. 1.—The radio index of different types of galaxy. 
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similar to those found in Sb and Sc spirals, otherwise we would expect their 
radio indices to be much lower than the radio indices of the Clouds. It follows 
that any corona which may exist in an irregular must emit less than the main body 
of the nebula. 

The above conclusion is supported by a comparison of the emission from the 
disk of a spiral and from the body of an irregular. Thus the radio index of an 
average Sb or Sc spiral, if we exclude about go per cent of the emission which is due 
to the corona, probably lies between + 3 and +4 and is therefore comparable with 
the value of +3 shown in Fig. 1 for irregulars. ‘This suggests that the radio 
emissivity of the disk of a late-type spiral is comparable with that of the body of an 
irregular galaxy and that the emission processes are similar in the two cases. 

The preseni limits to the radio index of early-type galaxies do not permit any 
positive conclusions to be drawn about their radio emission; nevertheless it 
appears likely that they are relatively weak radio emitters in comparison with Sb 
and Sc spirals. 


4. Acknowledgments.—We thank the Director of the Nuffield Radio Astronomy 
Laboratories for making available the facilities for this investigation. 
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NOTE ON 


THE GRAVITATIONAL INSTABILITY OF AN INFINITE 
HOMOGENEOUS ROTATING VISCOUS MEDIUM IN THE 
PRESENCE OF A MAGNETIC FIELD 


Since the publication of the above paper (M.N., 122, 455-459, 1961), 
Professor St. Piotrowski has kindly pointed out that equivalent but more 
extensive work on this problem has been carried out by his colleagues, 
A. G. Pacholezyk and J. S. Stodolkiewicz (Acta Astronomica, Warsaw, 10, 1, 1960). 
I should like to acknowledge the priority of their work, of which the Society and I 
were unaware at the time of publication. 

Similar work has also been carried out more recently by S. S. Kumar (Pub. 
Astron. Soc. fapan, 12, 552, 1960) and I am grateful to Dr Kumar for pointing 
out that the term 4y grad div u should be included on the right-hand side of the 
Navier-Stokes equation (1). Nevertheless, this does not affect the conclusion 
of my paper that the viscosity of the medium does not alter the Jeans criterion 
for instability. 
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